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CHAPTER I 


INTRODUCTION AND BACKGROUND 


A. Introduction 

The topic of high-intensity sound propagation has been mostly 
limited to sound waves in fluids and solids (see Beyer [lO] ). Very 
little has been done in the study of high-intensity sound propagation in 
bulk porous materials. Studies of intense sound interaction with per- 
forated sheets, thin sheets of porous materials, and Helmholtz resona- 
tors have been done (for example, see Refs. 37, 48, 55, 86). 

Recently the use of porous materials in jet engine inlets has 
increased. At high intensities it has been observed that the impedance 
of a porous material changes with sound level, i.e., the material be- 
haves nonlinearly. An understanding of sound reflection from and propa- 
gation in bulk porous materials has therefore become necessary. In this 
study many experiments have been performed on air saturated, fibrous, and 
expanded plastic porous materials. The materials have porosities in the 
range 0. 809 < P < 0. 985, where P is the porosity (the volume of air per 
total volume) . The emphasis of this study is on the experimental re- 
sults and explanation of these results. 

In this investigation a theoretical model for intense sound 
propagation in very porous, rigid, air saturated, fibrous bulk ma- 

* 

References are listed alphabetically by author and referred to by 
number at the end of this dissertation. 
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terials is developed. A separate empirical model is proposed to des- 
cribe finite amplitude losses (inexplicable losses in linear theory) in 
the porous material. Modeling of the reflection of high intensity 
sound waves from the surface of porous materials was also carried out. 
The reflection process is very complicated; both resistive and reactive 
components can depend upon the particle velocity amplitude of the in- 
coming wave. Experimental data obtained from measurements on a variety 
of bulk porous materials are used to test the validity and limits of the 
theoretical models. 

After this introductory section, some background theory is 
presented in the remainder of this chapter. First, a brief discussion 
of the theories on low intensity sound propagation in porous materials 
is used to put the theory of Chapter II in perspective. Second, a brief 
summary of high-intensity sound propagation in ordinary fluids is pre- 
sented. Third, past work on the reflection of high-intensity sound from 
surfaces and thin sheets is described. 

In Chapter II a high-intensity sound theory is presented. 

The one-dimensional mass, momentum, and internal energy equations for a 
nonlinearly behaving bulk porous material are derived. The resulting 
equations are solved by perturbation. A single boundary condition is 
used. The input signal is assumed to be a sinusoid which is distorted 
by a second harmonic component. A mathematical approximation to model 
the high-intensity impedance of semi-infinite materials is also pre- 
sented. An empirical model is adapted from the nonlinear acoustics of 
fluids to help describe the excess attenuation and approach to 
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saturation* of an intense sine wave in a porous material. The various 
theoretical results from Chapter II are compared to data in Chapter IV. 

In Chapter III the devices and methods used in the experi- 
ments are described. The devices are used to determine the porosity, 
material structure, flow resistivity, acoustic impedance, and acoustic 
propagation parameters (attenuation and phase speed) . Data are pre- 
sented in Chapter IV. 

In Chapter IV the experimentally determined data are compared 
to the theoretical results derived in Chapter II. The results of the 
porosity and dc flow resistivity measurements are presented and then 
used in the theory to predict the acoustical properties of the various 
materials. In the linear region the propagation parameter test results 
agreed with theoretical predictions in some cases and disagreed in other 
cases. The perturbation theory was found to be a poor predictor of the 
data at high intensities, porosities, and nonlinearities. As the 
porosity was reduced the agreement between measurements and predictions 
was better at higher intensities. The excess attenuation model predicts 
the excess attenuation of the fundamental over a large range of sound 
intensities, material nonlinearities, and porosities. Results of small- 
signal impedance measurements on finite and effectively semi-infinite 

* 

Saturation occurs when the finite amplitude losses become so large that 
no matter how much energy enters a porous material, only a specific 
amount of acoustic energy (saturation level) will arrive at some loca- 
tion within the material. The saturation level depends upon the dis- 
tance the wave travels, nonlinearity, and small-signal attenuation [8 l] 
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materials are presented and, in most cases, agree with theoretical pre- 
dictions. Results of high-intensity impedance measurements for effec- 
tively semi-infinite materials are presented and compared to the 
mathematical approximation presented in Chapter II. The advantages 
and problems of each model are discussed. 

In Chapter V the investigation of acoustic waves in porous 
materials is summarized, conclusions are discussed, and proposals for 
future work are presented. The appendices include a theoretical analy- 
sis of heat transfer in fibrous porous materials, assembly drawings of 
some of the devices discussed in Chapter III, and computer programs. 
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B. Linear Theory 


A review of the various theories of low-intensity propagation 
in and impedance of porous materials is presented in this section. 

The following discussion is intended to help orient the reader with 
respect to the vast literature on porous materials and the present 
study. The intent is not to cover all the literature but to give a 
general view. High intensity effects are postponed to Section C of 
this chapter. The linear theories are arranged in four categories based 
on the model used. After discussion, each theory is related to the pres- 
ent analysis. The four categories are the scattering model, the capil- 
lary tube-fiber motion models, the rigid frame model, and the lumped 
element model. 

1. Scattering Model 

Sound wave scattering from the fibers of the porous material is 
an approach that has been employed with varying success. The approach 
has the advantage that no empirical constants have to be determined 
from measured data. One minor drawback is that the theory is developed 
only for very porous fibrous materials. 

In 1910 Sewell [72] applied this method to fibrous materials 
containing a viscous gas. At that time there were no experimental re- 
sults for comparisons. In 1970 Kozhin [42-44] reworked Sewell’s theory, 
but made no comparison with available experiments. These theories yield 
predicted attenuations that greatly exceed the attenuations measured in 
this study. 
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In 1970 Attenborough and Walker [3] published a scattering 
theory for fibrous materials. They show accurate prediction of the ab- 
sorption coefficient at frequencies above 500 Hz, accurate prediction of 
the phase speed, and predicted attenuation exceeding measured attenua- 
tion at all frequencies. The predicted acoustic impedance is high. The 
theory includes both rigid and flexible porous materials. A more de- 
tailed discussion is given in Ref. 77. 

In 1976 Mechel [53,54] derived a long and complicated scat- 
tering theory in which he includes the scattering of viscosity, thermal, 
and density waves. He concluded that the theory was not as good a pre- 
dictor as some of the empirically related theories, such as Delany and 
Bazley's [25] analysis of Zwikker and Kosten's model [88] (discussed 
later in this section) . 

Scattering theory cannot be used to predict all the acoustical 
properties of a fibrous porous material, but in many cases the theory is 
adequate. In the present study, which is of high-intensity sound, this 
theory is not applicable because of the assumption that particle dis- 
placement is small in the porous material. It would be very difficult 
to generalize scattering theory to make it apply to high-intensity waves. 

2. Capillary Tube-Fiber Motion Models 

Two types of models are discussed in this section. They are, 
in most cases, intertwined with each other to such a degree that they 
cannot be separated. Most of these models include material motion in 
the acoustical analysis. The major difference between the models is 
the manner in which the viscous and thermal properties of the material 
are determined. In one case, the capillary tube model, the viscous and 
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thermal properties are defined by calculating the acoustical effects of 
a large number of adjacent capillary tubes. In the other case, the flow 
resistance model, the dc flow resistance of the material is used to de- 
termine the viscous properties. In both cases, one or more empirical 
factors must be included to match the theory to data. The following 
review is presented as historically as possible. 

The first attempt at theoretical analysis of the acoustical 
properties of bulk porous materials was done by Rayleigh in 1883 [66], 

He expanded his work in the second volume of his book [67]. Rayleigh 
assumed that the porous material could be modeled as an array of packed 
capillary tubes with the sound traveling axially along the tubes. 

Porous material structures are actually much different from a simple 
model of packed capillary tubes. He based his work on Kirchhoff’s [ 4 l] 
theory of viscosity and heat conduction effects on sound waves traveling 
in circular tubes made of perfectly rigid and heat conducting walls. 

Many others [6,8,11,25,73,74,88] have used Rayleigh’s approach as a 
basis for analysis and have used complicated schemes to predict experi- 
mental data. Some of these theories are discussed here. 

Adaptation of Rayleigh’s theory (see, for example, Zwikker and 
Kosten [88], Beranek [6,8], and Bies [ll]) requires that an empirical 
parameter called the structure factor be determined. The definition of 
the structure factor depends on the author. The most common definition 
of the structure factor is that it is a correction for the tortuosity 
the sound wave encounters as it travels in the material [6,8,11,88]. 
Another definition of the structure factor is that it is a correction 
for the effective acoustic air density caused by motion of the frame 
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[ 1 1 » 88] . Although the theories discussed in this section can be tied, 
in special cases, to the theory used in the present study, the structure 
factor is not used here. 

Zwikker and Kosten [88] have presented the most comDrehensive 
study of the acoustics of porous materials. They used Rayleigh's basic 
theory [67] as a starting place and introduced the structure factor 
and frame motion. In turn, their theory has become the starting point 
for the rest of the theories discussed in this sub-section. 

Beranek [8] started with the equations that Zwikker and Kosten 
[ 88 ] derived in their analysis. Except for the structure factor, 

Beranek T s analysis reduces to the rigid material model described in the 
next sub-section. Beranek measured the propagation parameters and im- 
pedance and, through his analysis, presented results illustrating the 
volume coefficient of elasticity of air in the porous material (heat 
transfer effects) and the effects of having the sample vibrate in its 
holder. A theoretical analysis of heat transfer effects is discussed 
in Appendix A of this study. Beranek f s measured data show a much smaller 
transition region between the two heat transfer states than was deter- 
mined here or by others [11,34,88]. In determining the sample vibra- 
tions in its holder, Beranek addressed a problem important to avoid in 
the experiments. If the sample holder holds the sample too tightly, or 
too loosely, the sample will tend to resonate at low frequencies and 
cause the experimental results to be erroneous relative to results de- 
termined from properly held materials. The resonance problem was con- 
sidered in the present study and the experimental results indicate that 
the problem was avoided. Beranek also discussed coupling of acoustic 
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and solid waves, but concluded that, at higher frequencies, the ma- 
terial becomes decoupled from the acoustic waves. 

Delany and Bazley [25] carried out an extensive empirical study 

of the acoustical properties of fibrous porous materials. They norma- 

* 

lized their data to the dimensional parameter f/o, where f is the 
frequency in Hertz and a is the flow resistivity in MKS Rayls/meter. 

The empirical relations they determined are useful over the specified 
range of 0.01<f /o<l .0. The relations are used to calculate the attenua- 
tion, phase speed, and impedance. Although not directly applicable to 
the present study, the equations are useful in checking measurements 
made at low intensities. 

Bies [ll] reviewed some of the above theories. He also dis- 
cussed measurement of flow resistance and porosity. Bies gave an ap- 
proximate relation for the structure factor. The discussion was not 
limited to fibrous materials. He presented data on cloths, fibrous 
materials, and fiberboard materials. 

Most of the above theories account for fiber motion (other 

than sample resonance). Lambert’s [49] analysis of Zwikker and Kosten’s 

r ® 

theory L 88 J indicates that above 130 Hz the air and frame of Kevlar 29 

aramid (du Pont de Nemours Company) are "decoupled".* Because f > 100 Hz 
in the present study, we are not concerned with frame motion. Support 

All the symbols are listed together at the beginning of this disserta- 
tion. 

^Kevlar®29 is a plastic fibrous material that is used extensively in the 
present study. 
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for Lambert's conclusion is given in Chapter IV. Lambert also deter- 
mined that the decoupling frequency for Scottfoara (Scott Paper Company) 

* 

was at 85 Hz. Scottfoam is a fully reticulated, expanded polyurethane 
foam and is similar to the Scottfelt (Scott Paper Company) used in the 

3 

present study. Although Lambert used the material density p g = 600 kg/m 
when, in actuality, the density of Scottfoam is p g =1153±32 kg/m [89], 
the error leads to an error of only 4 Hz in the decoupling frequency. 

We conclude from Lambert's analysis that the rigid frame model is appro- 
priate for use in this study. 

3. Rigid Frame Model 

In this section the low-intensity theory that is used as a 
basis and reference for the analysis in Chapter II is presented. This 
theory is based on the assumptions that (1) the dc flow resistivity can 
be used directly in the momentum equation and (2) the material structure 
is rigid. The first to use these assumptions were Kuhl and Meyer [45]. 
By using the flow resistivity a, their paper (1932) represents a simpli- 
fication of Rayleigh's more complicated treatment. Other authors [23, 
30,32,33,68] have used or rediscovered K\ihl and Meyer's approach as a 
starting place and expanded on it. 

Kuhl and Meyer used equations equivalent to the linearized 
continuity equation 

Sp >t +^=° * (I - 1} 

In the process of making a foamed plastic, thin plastic membrances form 
between the air cells of the foamed plastic. In a fully reticulated 
foam the membranes have all been removed. 
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the linearized momentum equation 


P u . 


P 



(1-2) 


2 . . 
and the equation of state P=c p to solve for the impedance, attenuation, 

and phase speed, where u is the acoustic particle velocity in the x di- 
rection, p is the acoustic pressure, Sp is the acoustic density, p Q is 
the ambient density, c is the sound speed, P is the porosity, x is the 
distance, t is the time, and the comma denotes differentiation of the 
dependent variables with respect to x or t (the independent variables). 
In their theoretical analysis Kuhl and Meyer assumed the compressibility 
of air in a porous material to be variable. They asstimed the adia- 
batic sound speed c 0 = A'I > 0 /p o ? where y=Cp/C v Is the ratio of specific 
heats. They reasoned that the compressibility, and thus the density, of 

the air could be determined by experiment. In their calculations they 

2 

used the standard compressibility of air Q q C q . In Appendix A we assume 
a constant density P Q * We show that, because of heat transfer effects, 
the speed of sound is a function of both porosity and frequency. Be- 
cause the relation derived in Appendix A is frequency dependent, we use 

the isothermal sound speed b = /p 7p . 

o 00 

The above equations may also be combined to obtain a wave 
equation [23] 


u 

,xx 
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b 2 P 

00 




(1-3) 
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The use of P as a divisor of u depends upon whether the reference u is 
in the material or in the open air (no material). Here, as P decreases 
the particle velocity in the material increases. The importance of 
using u/P as the effective particle velocity is discussed in Chapter II. 
Note that although Eq. 1-3 seems to be independent of P, the flow re- 
sistivity a is actually a function of P. 

To solve Eq. 1-3, one can assume that a time harmonic wave 


j (wt-rx) 

u=e 


where T^g-ja, w=2Trf, and j=/^T, propagates in a porous material. The 
wave encounters the attenuation 


a 



(1-4) 


and the wave number 



The phase speed c =u)/g is 

PH 

7i b 


o 



(1-5) 


( 1 - 6 ) 


Plots of attenuation and of phase speed versus frequency for several 
commonly found resistivity values are shown in Pigs. 1-1 and 1-2, 
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FIGURE 1-1 

ATTENUATION VERSUS FREQUENCY FOR THREE FLOW 
RESISTIVITIES OF A POROUS MATERIAL. 
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FIGURE 1-2 

PHASE SPEED VERSUS FREQUENCY FOR THREE FLOW 
RESISTIVITIES OF A POROUS MATERIAL. 




respectively. Both the attenuation and phase speed can change dramati- 
cally over a relatively small frequency range. These changes are im- 
portant in the propagation of intense sound in porous materials. In 

porous materials intense sounds do not behave in the same manner as in 

2 

open air, which has a low, uniformly varying attenuation (a«u> ) and only 
minor (usually negligible) dispersion. 

Kuhl and Meyer showed that the specific acoustic impedance in- 
ternal to an essentially semi-infinite (aL>>l, where L is the sample 
length) porous material is 


w = W R + j w i 


where the resistive part is 



and the reactive part is 


W = 
I 



(1-7) 


( 1 - 8 ) 


(1-9) 


We see that Eqs. 1-4 and 1-5 are of equivalent forms to Eqs. 1-9 and 
1-8, respectively. We conclude that, for the small signal case, 
measuring the impedance when aL>>l is equivalent to measuring a and 0. 
This comparison proves to be a good check on both the theory and data. 

The real and imaginary parts of the impedance are plotted 
versus frequency for three flow resistivities in Fig. 1—3. 
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FIGURE 1-3 

ACOUSTIC IMPEDANCE VERSUS FREQUENCY FOR THREE 
FLOW RESISTIVITIES OF A POROUS MATERIAL. 



are 


The general propagation parameter T and impedance W equations 
related to each other by [49,88] 

W(L) = (W +Jfty coth(a + jB)L . (1-10) 


Equation 1-10 is for finite length (aL < 1) materials mounted on a hard 
surface. Equation 1-10 is appropriate for most low-intensity impedance 
measurements and is used in the analysis presented in Chapter IV. 

When a sound wave traveling in air encounters a porous material 
some of the energy of the incident wave is reflected from the material 
surface and some is transmitted into the material. The Impedance dif- 
ference between the air p c and the material W(L) causes the separation 

o o 

of the incident wave into two parts. The ratio of the transmitted (ab- 
sorbed) energy to the incident energy is the absorption coefficient. 


a 

n 


W(L) - p c 
o o 

W(L) + p c 
o o 


( 1 - 11 ) 


As the transmitted wave travels in the material it is attenuated (Eq. I— 
4) by the action of viscous effects. In the present study the term ab- 
sorption is used to describe a reflection/transmission process, whereas, 
the term attenuation is used to describe the decay of a traveling sound 
wave. Absorption coefficients (aL>>l) are plotted versus frequency for 
three flow resistivities in Fig. 1-4. 

In Figs. 1-1 through 1-4 we find that the value of the flow 
resistivity affects the acoustical properties of porous materials. Thus, 
the importance of determining flow resistivity cannot be overstated. The 
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FIGURE 1-4 

ABSORPTION COEFFICIENT VERSUS FREQUENCY FOR THREE 
FLOW RESISTIVITIES OF A POROUS MATERIAL. 



flow resistivity a depends on the material structure. Flow resistivity 
may be determined either through acoustical testing or through dc 
(steady state) flow resistance testing. Both methods are used in the 
present study. At this time, there are no theoretical relations that 
can be used to predict a. Several empirical relations [24,33,34,60,78] 
have been formulated to predict the flow resistivity of fibrous materials. 
Only two relations [33,34] will be discussed in the present study. 

By using hydrodynamic analysis, Hersh and Walker [35] deter- 
mined a relation between viscosity, porosity, fiber diameter d, and flow 
resistivity inside a fibrous porous material with flow perpendicular to 
the fibers, 

a - 



where g is an empirical constant. Equation 1-12 is used in Chapter IV 
to calculate the dc flow resistivity for unmeasured values of a. 

Hersh and Walker determined g=0.059 for batted Kevlar 29. They com- 
pared data and theory at low frequencies. 

For their data-theory comparisons Hersh and Walker derived a 

low frequency approximation to Eq. 1-3 by ignoring the u term. For 

, tt 

a sinusoidal wave they found that the attenuation a and wave number £ 
were equal to 



(1-13) 
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For porosities in the range of 0.9<P<0.97, the low frequency approxi- 
mation is valid for f<2 kHz. For a full frequency range prediction the 
constant g must be redefined. 

Hersh and Walker [33] showed that a is related to d, y, and P 


by 



(1-14) 


where h(P) is a structurally dependent function of P. For their specific 
case the relationship is given by Eq. 1-12. In dimensionless quantities 
Eq. 1-14 is 


a = y 2 h(P) = h(P)/Re , (1-15) 

p wd 
o 

where Re is an acoustic Reynolds number [5,33]. The value of the acous- 
tic Reynolds number is important in our analysis. When Re>>l the acous- 
tic boundary layer is small relative to the fiber size, viscous effects 
are unimportant, and Eq. 1-12 is not valid. When Re<<l the boundary 
layer is large, viscous effects are important, and Eq. 1-12 is valid. 

For the materials and frequency range used in the present study, 

-3 

10 <Re<23. Most of the experiments were done such that Re< 1 . 

Flow resistance data for many different porous materials have 
been presented in several studies. Nichols [60] dealt solely with the 
linear flow resistance of fibrous porous materials. He also determined 
an empirical relation for fibrous materials. Brown and Bolt [19] pre- 
sented data on both linearly and nonlinearly behaving materials and were 
the first to plot flow resistance data versus particle velocity on a log- 
log plot. Plotting the data in this manner illustrates the relative 
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nonlinearity of the dc flow resistivity (see Chapter IV) . Bies and 
Hansen [12] present linear flow resistance data for many porous materials. 

Only viscous effects are considered in the above theory. The 
frequency dependent effects of heat transfer on the compressibility of 
the gas in the porous material are ignored. Heat transfer effects 
cause the compressibility of the air to be neither adiabatic nor iso- 
thermal and, for comparison with data and for simplicity in the calcula- 
tions, the limiting values of the adiabatic and isothermal sound speeds 
are each used in the low-intensity sound calculations. In the high- 
intensity sound calculations, because of increased heat conduction, only 
the isothermal sound speed is used. The theory that has been summarized 
up to this point will be used as a basis for comparison between the first 
order theory developed in Chapter II and the small signal experiment re- 
sults presented in Chapter IV. 

In 1980 Hersh and Walker [33] extended their 1979 theory [34] 
to include heat transfer effects. They used results of experiments on 
dc flow through fiber bundles to determine relations to fit equations 
for heat transfer and viscosity effects. They obtained a one-dimensional 
wave equation that combines both heat transfer and viscosity effects. 

There are three empirical constants to fit to the data; one constant 
deals with heat transfer and two constants deal with viscosity. The 
values of the viscosity constants depend on whether the sound travels 
perpendicularly or axially with respect to the fiber orientation. The 
wave equation is 
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where K is a heat transfer parameter and a is a viscous drag parameter; 

n 

both parameters are functions of the porosity, fiber size, and their 
respective fitted constants. The viscous drag parameter presented by 
Hersh and Walker [34] is 


4y 0 lii [1/ D + V D ], 
d 2 n n p p 


where 

D = 16/TTp [1 + 14.75(1 -P) 3 ] , 
n 

D = 3.94(1 -P)°’ 413 [l + 27(1 -P) 3 ] 
P 


(1-17) 

(I-18a) 

(I-18b) 


where l / and 1/ are the empirical constants and the subscripts denote 
n p 

flow parallel (p) to the fibers and flow normal (n) to the fibers. Hersh 
and Walker fit Eq. I-18a to dc flow resistance data taken on many dif- 
ferent materials by Davies [24], Davies found very little scatter of 
measured data from his own theoretical predictions. Hersh and Walker 
fit Eq . I-18b to dc flow resistance data taken on a variety of compact 
bundles of parallel fibers by Sullivan [78]. Sullivan also found little 
deviation of measured data from his own theoretical predictions. Even 
though Davies and Sullivan found little scattering in their data, the 


data fitting constants l/^ and V in Eq. 1-17 


should be determined for 


each material to be used. It appears that, since there are two constants. 
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the equation may be used for a material made of randomly oriented fibers 

[34], As shown by Hersh and Walker and in Chapter IV, the definition of 

a can be used to define the flow resistivity without including heat 
H 

transfer effects. 

The heat transfer parameter presented by Hersh and Walker is 

K * ^■6YU(l-P) V2 [lO- 94Il^Pjil _ (I _ 19) 

Prd PT 

n 

where Pr=yC /< is the Prandti number, k is the thermal conductivity of 
P 

the medium, and T is an empirical constant. The basis for this 
n 

equation was taken from work done by Masliyah [55] on dc energy and 
momentum transfer from cylinders oriented normally to the fluid flow. 

The heat transfer effects on the sound speed in Eq . 1-19 are similar to 
those determined in Appendix A. 

The material Hersh and Walker used in their tests is Kevlar 29. 
As noted earlier, this material is used. in many of the experiments in 
the present study. In addition, other materials are used in the present 
study and the test results are discussed in Chapter IV. 

4. Lumped Element Model 

The last linear theory to be considered here was published by 
Zarek in 1978 [82], He considered sound absorption by flexible poly- 
urethane foams. He considered materials with and without an Impermeable 
membrane facing. He approached the theory from a lumped parameter model 
for the gas and material and used the results to determine a Lagrangian 
that includes Rayleigh dissipation. He then applied Hamilton's 
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principle and determined equations of motion for the coupled 
air-material system. He stated that the parameters only include the 
small scale effects of the porous medium; all nonlinearities of the air 
and material are ignored. His theoretical predictions agree well with 
his experimental results. Zarek also compared his theoretical pre- 
dictions to some of Beranek's [7] data and found fairly good agreement 
between the two. Zarek T s theory is not useful in the present study because 
he considered material motion and he ignored all nonlinearities. 

The theories discussed in this section all pertain to the 
acoustical properties of the porous materials at low intensities. Many 
of the theories include frame motion of the material, but Lambert T s 
analysis shows that, for the present study, the frame motion is negli- 
gible. Neglecting the frame motion greatly simplifies the analysis. 

In Chapter II we use the conservation equations and show which assump- 
tions lead to the rigid frame theory. 
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C. Nonlinear Theory 


This section contains a presentation of two topics. The first 
topic is a brief review of how an intense sound propagates in a fluid 
and how absorption and dispersion affect shock formation and attenuation 
The second topic is a brief presentation of the theoretical and experi- 
mental studies of intense sound reflected from thin porous sheets. 

1. Intense Sound in a Fluid 

The propagation of an intense sound in a fluid is fairly well 
understood. The brief theoretical treatment presented here is based on 
extensive work by others and is presented as background to the nonlinear 
theory used in Chapter II. 

For any point on the waveform the propagation speed will vary 


as [lO] 


dx 

dt 


u=const 


c + u 


( 1 - 20 ) 


Equation 1-20 adds the effect of convection u to the sound speed 


c = 


/ dP 


/dp 

s 


( 1 - 21 ) 


where P is the total pressure, s is the entropy, and p is the total 
density. The pressure-density relation is nonlinear, c depends on tem- 
perature, and for a simple wave 


c 


c +^ u 
o 2 


( 1 - 22 ) 



The propagation speed becomes 

= c +8’u , (1-23) 

_ o 
u=cons t 

where g'=(y+l)/2 is the nonlinearity parameter of air. In the small 
signal case, where u<< c q , 

dx 

dt C o 

u=const 

In an isothermal situation, such as in a porous material, y=l and 6 ,= 1 
and the nonlinearity of the air is reduced. In a porous material the 
nonlinear effects are more complicated than suggested here. 

Since the propagation speed varies with location on the wave- 
form the wave peaks will travel the fastest and the wave troughs the 
slowest. The thermal and convective effects are cumulative, the wave- 
form distorts as it travels, and, when u is large enough, the compres- 
sional phases steepen and form shocks. Distortion of the wave is accom- 
plished by the generation of higher frequency harmonic components. All 
fluids dissipate energy, the high frequency components are attenuated at 
a faster rate than the fundamental component, and, thus, dissipation 
limits the distortion of the waveform. By shifting energy to the higher 
harmonic components, the fundamental loses energy more rapidly than at 
low intensities. The increased loss is called excess attenuation and we 
find that the excess attenuation increases with source intensity. Above 
a certain source level we find that, no matter how much we increase the 


dx 

dt 
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energy input, the excess attenuation limits the measured level to a 
constant. We call this constant level the saturat ion level [A , 57 , 76 , 8l] . 
The farther a measurement point is from the source the lower the satura- 
tion level is at that location. Before the dissipative effects are dis- 
cussed, the progressive distortion of the wave is discussed. 

In studying shock formation and dissipation in a lossless 
fluid, the Fubini solution [16] 


f-E 


2 * 

— J (n a )sin n(cot-kx) 
* n 
n a 


(1-24) 


applies to the region before shocks are formed (ct*< 1) when the boundary 
condition at x=0 is 

u(0,t) = u^sin cot , (1-25) 

th 

where u is the initial particle velocity amplitude, J (•) is the n 
o n 

order Bessel function, a*=x/x, x=l/3 ? ek is the shock formation distance, 

e=u /c , k=aj/c , and n denotes the n ^ harmonic, 
o o o 

In the region (cr*>3) where shocks are well formed and have be- 
gun to decay, the Pay solution [16] 

— = ^ ^ ^ * sin n(oot-kx) (1-26) 

U o sinh[n(l+a )/g] 

is applicable, where G=B T ek/ct, which is sometimes called the Gol'dberg 
number [lO], In the limit, as G-*°°, Eq. 1-26 represents a sawtooth wave. 
In this limit the ordinary dissipation is small relative to the rate at 
which energy is pumped into the higher harmonic components and dissi- 
pated In the shocks. Blackstock [l6] devised a scheme to combine the two 
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solutions, Eqs. 1-21 and 1-23 (G=°°) , into a third solution such that 
the whole region from the boundary to infinity may be explored. 

Since G^°°, a perfect sawtooth wave never forms. The Gol T dberg 
number indicates the relative importance between distortion and ordinary 
dissipation. For a very strong wave, such that B'ek>>a, the shock forms 
quickly, approximating a sawtooth, and takes some distance to dissipate. 
For weak waves or large viscous dissipation, such that B'ek<<ct, a shock 
wave will not form. 

In 1977 Bjcirn^ [l4] compared dissipation effects to nonlinear- 
ity effects by using both the Gol'dberg number and the Keck and Beyer 
[40] perturbation solution of a nonlinear wave equation to study in- 
tense sound propagation in a viscous fluid saturated, rigid sediment. 
Bj$Srn?5 experimentally determined the parameter of acoustic nonlinearity 
in fluid-saturated sediments by using Beyer ? s analysis [9], Bjcirn^ 
showed in his analysis that the amplitude of the second harmonic com- 
ponent could not exceed 1.3% of the fundamental components because the 
attenuation effects are much larger than the nonlinearity effects. 

Bj^rn^ indicated that, for most high-amplitude sounds, the sound goes 
directly from finite-amplitude distortion into the "old age" (the wave- 
form is nearly sinusoidal) propagation region. He noted that a shock 
wave propagating in a sediment decays rapidly and he concluded that the 
nonlinear effects are negligible in viscous fluid saturated sediments. 

Determining the nonlinear behavior of a material is important. 
If the material behaves very nonlinearly, relative to attenuation, then 
harmonics are generated very rapidly. In the present study the dc flow 
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resistivity characteristics are used to define both the viscous dissi- 
pation and nonlinear properties of the material. Carman [ 2 1 ] stated 
that for dc gas flow through porous materials the transition to non- 
linear behavior is related to the onset of turbulence and the modified 
Reynolds number 

R - A , (1-27) 

m vo 

where v is the kinematic viscosity and S=4(l-P)/d is the surface area/ 

unit volume. Carman stated that when R > 1 the flow resistivity behaves 

m 

nonlinearly, i.e., the flow resistivity depends on the particle velocity. 
Later in this study we find that the manner in which the nonlinearity is 
caused, i.e. , turbulence, may be important to how an intense sinusoid 
loses energy. We find that the complexities of excess attenuation of a 
sinusoid are not as easily defined in a porous material as in an open 
fluid. In addition, how the nonlinearity is defined and used becomes 
extremely important to how well the effects of nonlinearity are pre- 
dicted. 

Blackstock [ 1 7 ] presented a perturbation solution for a 
Burgers equation governing the propagation of sound in an absorbing, 
dispersive fluid. He showed that the dispersive characteristic affects 
the level of the generated second harmonic component. This effect has 
importance in a porous material. Equation 1-7 shows that porous ma- 
terials are very dispersive. The harmonic components of intense sounds 
generated in a porous material may be greatly affected by the dispersion 
because energy added to the components will be added out of phase with 
the energy already present. 
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Zorumski [84] presented conservation equations for thin sheets 
of porous materiais with nonlinear dc flow resistance. He used these 
equations to determine the scattering of sound from porous elliptic 
cylinders. He did not attempt to solve the conservation equations for 
the propagation of intense sound in a bulk porous material. His re- 
sults are described in the next section. 

In this section a very brief review on intense sound prop- 
agation through lossless and dissipative fluids has been presented. 

These ideas are useful in the analysis and comparisons made in Chapter IV. 

2. Reflection of Intense Sound 

For all practical purposes, a sound propagating in a porous 
material must first enter the material from an adjacent medium. A wave 
incident on the surface of a bulk material is reflected from as well as 
transmitted into the material. In order for us to understand the losses 
associated with the transition from one medium to another, we must deter- 
mine the impedance of the material. Most materials are acoustically 
fi n ite (otL< 1) and, as shown in Eq. 1-10, the impedance, the propagation 
parameters, and the material termination impedance (at x=L) all influence 
the measured impedance. At high intensities the analysis becomes dif- 
ficult because both the impedance and propagation parameters depend on 
intensity . We have just seen that as an intense wave propagates, energy 
is transferred to the higher harmonic components. The harmonic components 
interact with each other and superposition is no longer valid. When 
an intense wave is used to measure reflection (or absorption) properties 
of a porous material, the sound level, instead of just frequency, must 
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also be specified. Since superposition is not valid, the incident and 
reflected waves may influence each other. For short interaction regions 
this influence may be small. For large interaction regions, such as in 
long pulse trains or standing waves, the mutual influence must be ac- 
counted for. 

This section contains a review of past studies of the reflec- 
tion of intense sound from different surfaces and the determination of 
the impedance of sheets of porous materials. Both theoretical and ex- 
perimental studies are examined. 

In 1960 Blackstock [ 1 5] showed that intense sounds reflected 
from surfaces do not always follow the commonly used small signal laws. 
Blackstock considered a pressure release surface, an infinite impedance 
surface (a hard wall), and a thin resistive sheet. Regardless of in- 
tensity, reflections from a pressure release surface are accompanied by 
a doubling of the particle velocity at the interface. In addition, 
Blackstock found that at a hard wall the variational sound speed c is 
the variable that doubles throughout the intensity range. The pressure 
doubling law for reflection from a hard wall is only an approximation. 
However, at a sound pressure level of 174 dB re/20 yPa the deviation 
from pressure doubling is only 6%. Above 174 dB pressure amplification 
at the wall increases and a more exact relation must be used. 

In the case that most interests us here, the thin resistive 
sheet, Blackstock found a relation for the pressure reflection coeffi- 
cient 
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where P^ Is the reflected pressure, P^ Is the incident pressure, 
4(z)=Re(z)-p o C o , z is the specific normal acoustic impedance, lm(z)=0, 
e=u /c is the acoustic Mach number of the initial wave, and 6 is a 
specific location on the wave that is followed through the course of 
travel of the wave. 

Unfortunately, porous materials do not meet the condition that 
lm(z)=0. However, the case of purely resistive porous materials is in- 
structive. From the analysis results the acoustical characteristics of 
the real material may be inferred. Equation 1-28 may be used to follow 
only one point on the waveform at a time. If the first peak of the sine 
wave (compression) is followed as it leaves the piston, sinu)<t> = l . 0. In 
terms of the real part of the acoustic impedance, the absorption coeffi- 
cient is 
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Equation 1-29 can be used to calculate the absorption coefficient that 
the peak of an initially sinusoidal wave encounters. 


32 



If the minimum of the sine wave (rarefaction) is followed as 

it leaves the piston, sinw<f)=-l . 0 and the signs of the two e terms in 

Eq. 1-29 become negative. This change of sign indicates that the phase 

of the wave is important to reflection from purely resistive porous 

materials. This statement of phase dependence should hold for any 

porous material. For example, when Re(z)/p c =3.0 and e=0.02 (160 dB 

o o 

re/20 yPa) a =0.763 for sin<f> = 1.0 and a =0.738 for sind>=-1.0. The dif- 
n n 

ference increases with increases in Re(z) and e. 

In the above analysis Re(z) is arbitrary and may be a function 

of the particle velocity. If the impedance is a function of the particle 

velocity and lm(z)^0, then a^ changes more dramatically than indicated 

by Eq . 1-29, because both Re(z) and Im(z) affect the value of a . 

n 

Zorumski and Parrot [83,86] and Zorumski [84,85] are the only 
researchers to present a combined theoretical-experimental analysis of 
thin nonlinear ly behaving porous materials. The most general presenta- 
tion is in Ref. 86. The mathematical and computational details are in 
Refs. 84 and 85, respectively. We summarize their general theory, ex- 
periments, and results below. 

In this theory Zorumski and Parrot present generalized, func- 
tional relations for the conservation of mass and momentum. They assume 
the material is thin enough that X >> L, where X is the wavelength, and 
that the velocity differential across the sample is negligible. They 
also assume negligible compressibility in the material. Thus, the con- 
servation of mass equation does not enter into their analysis. Their 
analysis is based on use of the momentum equation. They assume that the 
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particle velocity dependent dc flow resistance can be used to describe 
the acoustic flow resistivity R[v(t)], where V(t) is the acoustic 
particle velocity. In addition, they assume the acoustic reactance 
X[v(t)] is a function of particle velocity. From their version of the 
momentum equation they obtain a temporal impedance operator Z^_ such that 

Z t V(t) = |r[v(0] + x[V(t)] ^} v (t) • (1-30) 

The temporal impedance operator is evaluated from their experimental 
data. They claim the theory can be used to evaluate the impedance en- 
countered by any distorted waveform. Their analysis is done for dis- 
torted sinusoids. As discussed below, if there is no nonlinear inter- 
action between the harmonic components nor between the incident and re- 
flected waves, then their claim appears to be valid. 

We now discuss Zorumski and Parrot’s [83,86] experimental in- 
vestigation. They used a standing wave tube with a microphone flush 
mounted with the surface of the rigid termination of the tube. They 
placed the material at one-fourth wavelength from the rigid termination. A 
second microphone was placed outside the material surface that was away 
from the termination. Tests were done at 0.5, 1.0, 1.25, 2.0, and 
4.0 kHz and at sound levels from 120 to 160 dB. They noted, that, be- 
cause of the material and second microphone placement, all even harmonic 
components of the fundamental were ignored in the experiment and, thus, 
in the theory. Ignoring the even harmonic components might lead to 
problems because the components may interact with the waves everywhere 
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except for in the material. As seen in Eq. 1-24 harmonic component gen- 
eration from a sinusoid starts by generating all harmonics, and, mainly, 
the second harmonic. If all the even harmonics are ignored, then much 
of the energy lost by the fundamental is being ignored and measurement 
error is being introduced. 

Zorumski and Parrot presented several interesting results. 
First, they showed that the acoustic flow resistivity is independent of 
frequency and closely approximates the dc flow resistance at all par- 
ticle velocities. Second, they showed the acoustic reactance to be a 
function of both frequency and particle velocity. The theoretical- 
experimental results are found to be consistent and in fairly good 
agreement over large frequency and sound level ranges. 

The above analyses presented by Blackstock and Zorumski and 
Parrot show that the analytical and experimental determination of non- 
linear porous material impedance is extremely complex. This complexity 
has led many researchers to experimentally determine the impedance and 
absorption of various materials. Shock tube experiment results [27,35, 
79] will not be discussed. 

In 1970 Powell and Van Houten [64] used band-limited tone 
bursts at frequencies between 500 Hz and 10 kHz to study the absorption 
properties of porous material covered resonators. They measured the in- 
cident and reflected peak pressures at a single microphone position and 
discussed design considerations for the test procedures. The evaluation 
of their results is difficult because the exact length and diameter of 
their wave tube, or tubes, was not specified. It appears that the tube 
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length could have been greater than the shock formation distance at 
high frequencies and sound levels. They discuss design criteria for 
the optimum tube lengths and diameters, but do not state how, or if, the 
criteria were implemented. They analyze the waveform and spectral con- 
tent of what appears to be their input waveform to the acoustic driver. 
They never indicate what the incident and reflected acoustic waveforms 
look like. They did not indicate whether they saw shock formation or if 
the acoustic waveform was distorted by the reflecting surface. The tube 
length is important because at high intensity, if the tube length is too 
long, then the measured absorption coefficient would be incorrect. 
Measurement of the absorption coefficient depends on the attenuation of 
the traveling wave being constant with sound level change. A wave 
traveling in a tube longer than the shock formation distance x ex- 
periences different attenuations. The attenuation depends on where the 
wave is with respect to x, Before the wave reaches x, the shock is 
forming and the attenuation of the fundamental component is greater than 
tube wall attenuation. After x, the attenuation increases and then, 
when x>4x, asymptotically approaches tube wall attenuation [l6]. Some 
nonlinear attenuation effects can be accounted for in the calibration 
measurements, but, as indicated above, the error can grow rapidly and 
a short tube is necessary. 

In 1973 Melling [55] experimentally and theoretically deter- 
mined the impedance of perforated plates and perforated plate resonators 
at low and high intensities. He showed how both the real and imaginary 
parts of the impedance change with the intensity. He designed a high 
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intensity standing wave impedance tube which has a frequency range of 
300 to 3250 Hz and upper sound pressure limits of 144 to 159 dB [56], 
The system is large and also requires large transducers to generate the 
required steady state sound levels. Melling’s program is well thought 
out and could be useful for evaluating bulk porous absorber character- 
istics. Unfortunately, the cost of building a system of this type is 
large and this type system was not used in the present study. 

In 1980 Nakamura et al., [58] described the absorption of a 
small-signal, plane N wave by finite length porous materials. They 
showed that the individual Fourier components of the N wave interact 
individually with the porous material. If the material is long enough, 
one can determine and and, if it is short enough (aL<l), one can 
determine a and 8 . 

In a later paper Nakamura et al. [59] described the reflec- 
tion of a plane N wave from the end of an open pipe. They showed, as in 
their earlier paper, that the low level impedance is predicted by low 
level N wave tests, but the high intensity N wave reflections are 
drastically modified as the amplitude increases. They stated that their 
results show the total reflected energy between zero and 10 kHz to be 
approximately Invariant with respect to the N wave amplitude. In this 
latter paper, they used a much different system in the measurement of 
impedance. The large amplitude N wave does not propagate in a linear 
fashion and the frequency components interact. They used an algorithm 
similar to the algorithm developed by Pestorius [62] to mathematically 
propagate the measured incident wave to the pipe termination. They 
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modified Pestorius* algorithm to mathematically propagate the measured 
reflected wave backward from the microphone to the pipe termination. 
After they mathematically propagated both waves to the end of the pipe, 
they subtracted the Fourier transform of the reflected waveform from 
the Fourier transform of the incident waveform. They used the results 
of this subtraction to calculate the pressure reflection coefficient and 
the energy reflection coefficient of the open pipe. This method works 
because the reflected wave is inverted with respect to the incident 
wave and much of the high frequency energy of the incident wave is lost 
out the end of the pipe so that no shock exists in the reflected wave. 

In fact, by the time the reflected wave reaches the microphone, the wave 
does not develop a shock. 

If a porous material was placed in the system of Nakamura 
et al. [59], the reflected wave would (usually) not invert. A shock 

could reform in the reflected wave and the above method could not be 
used. Once a shock has formed in a wave the process cannot be mathe- 
matically reversed to obtain the original waveform because information 
has been irretrievably lost in the shock formation process. The method 
devised by Nakamura et al. Is limited because it may not be useful 

above certain sound levels and requires a computer for the calculations. 
For these reasons a more useful method was sought. 

In 1981 Kuntz et al. [47] described a simple impedance 
measurement method. This method is described in detail in Chapter III 
of this study and experimental data is presented in Chapter IV. 
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In this section several test methods and theoretical descrip- 
tions that show the change of impedance with intensity have been dis- 
cussed. The advantages and disadvantages of the different methods have 
been discussed and serve as background information for the measurement 
methods described in Chapter III. 
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D. Summary 

Several diverse topics have been considered in this chapter. 

In spite of their diversity, the topics are the basis for understanding 
the theory and results that follow in the remainder of this study. 

Many theories have been developed to explain low-intensity sound prop- 
agation in porous materials. In this study, because of the complexity 
of the nonlinear acoustics analysis, only the simplest approach is used, 
that of using the dc flow resistivity to define the viscous effects in 
a rigid porous material [30,32,33,45,68]. The more complicated ap- 
proaches that include heat transfer and material motion [8,11,34,49,67, 
88] could, conceivably, be used, but the mathematical complications are 
prohibitive. In most cases, the simple approach is seen to result in 
adequate prediction of the acoustical properties. 

From the nonlinear acoustic theory for fluids we have shown 
that, because of the high dissipation In a porous material, shock waves 
will probably not be formed in an intense sound wave. Perturbation of 
the conservation equations Is used in Chapter II to describe the prop- 
agation of intense waves in bulk porous materials. 

Propagation of a wave in the material is not the only aspect of 
sound interaction to be considered. The wave must enter the material 
and, consequently, a wave incident on the surface of a bulk porous ma- 
terial is reflected from as well as transmitted into the material. The 
mathematical prediction of the impedance that an incident, high-intensity 
sound wave encounters at the surface of a bulk porous material is not a 
trivial task. In fact, the task has never been accomplished. In this 


40 



study only a simple mathematical approximation to the impedance will be 
made . 

Because superposition is invalid for intense sound, the 
measurement of the impedance of a bulk porous material exposed to an in- 
tense sound is also not trivial. Many measurement methods have been 
developed to measure the impedance of both bulk and sheet materials. 
These methods all show that the acoustic impedance and absorption co- 
efficient of a porous material depend on the amplitude and frequency 
content of the incident sound wave. Another measurement method is pre- 
sented and used here because it appears to circumvent some of the pro- 
lems that the other methods do not address. 
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II. THEORY AND MODELING 


In Chapter I various theories for sound propagation in 
acoustic materials were described. The rationale for modeling porous 
materials as rigid materials was presented as a reasonable assumption. 

In Sections A through C of this chapter the mass, momentum, 
and energy conservation equations are derived for porous materials. The 
model is limited to rigid, isotropic, nonlinear, air saturated, bulk 
porous materials. 

In Section D perturbation is used to determine approximate 
solutions to the conservation equations. These solutions are used in an 
attempt to predict the manner in which intense sound propagates through 
a porous material. Because the measured boundary conditions are not 
ideal, an approximate boundary condition is considered. The boundary 
radiates both the first and second harmonic components into the 
material. The two components can be set to arbitrary amplitudes and 
relative phase. A second-order approximation is determined. In 
Chapter IV the solution is compared to the data. 

In addition to sound propagation in the porous material, the 
characteristic impedance of the material is important. In Chapter I we 
found this problem to be extremely complicated at high intensities. 

In Section E an approximation to the impedance relations is made to 
determine a useful intensity dependent impedance relation. The results 
only apply to sinusoidal waves normally incident on the surface of a 
semi-infinite porous material. 


42 



In Section F the approximation of Section E is extended to 
illustrate nonlinear propagation of the fundamental component. In 
addition, another model is proposed to illustrate the effects of 
saturation and excess attenuation on the fundamental component. 

In Section G a short summary of the chapter is presented. 
In Chapter IV the various models are compared to the data. 
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A. Continuity Equation Derivation 


The physical concept of how a fluid flows through a porous 
material is necessary for us to understand how the flow is distorted 
and impeded by the material. The continuity equation derivation is 
presented in simple physical arguments so that we may take certain 
assumptions for granted in future derivations. 

A simplified cross section drawing of a bulk porous material 
control volume is shown in Fig. 11-1, The structure is made of many 
parallel fibers. The control volume is fixed to the rigid structure. 
If no structure is present, the particle velocity is represented by u. 
When the fluid enters the porous material the velocity is increased by 
a factor of the porosity P. The porosity is the volume of the air in 
the control volume relative to the total volume (including fibers) of 
the control volume. We assume the material to be homogeneous and the 
porosity constant in each derivation. 

For the derivations we need to know the relation between po- 
rosity and the cross sectional area a fluid encounters in a material. 
In Fig. II-l a single "cell" has been drawn such that the four sides 
each bisect two fibers. Thus, a quarter of each fiber is enclosed by 
the cell. Assume each cell to be of length L along the fibers. The 
porosity of the cell is 


p = i - 


7T d 

4DHL 


where D and H are defined in Pig. II-l. 


(II-l) 
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FIGURE 11-1 

STATIONARY CONTROL VOLUME USED IN THE DETERMINATION OF THE 
CONSERVATION OF MASS RELATION FOR POROUS MATERIALS. 
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A gas flowing through the material will not have a constant 
speed. Thus, the average open area encountered by the gas flow may not 
be evaluated at any point in the medium. The open area, as a function 
of x, is integrated along the length D of the cell. In doing the in- 
tegration we find the same relation for the area as for the porosity, 

Eq . II-l. We define the average open area to be 

A p = PA , (II-2) 

where A is the cross sectional area of the control volume. 

A physical statement of the continuity equation is that 

the mass flux across the = the time rate of change of 
control volume surface mass in the control volume. 


For the one- dimensional system shown in Fig. II-l, the continuity 
equation is 

p £ Pa | — P ~ Pa 1 = Paax , (ii-3) 

p 1 x P 1 x+Ax 3 1 


where p is the total density, x is the position. Ax is the control 
volume length, and t is the time. A differential form of Eq. II- 3 is 
found by dividing Eq. II-3 by PAAx and taking the limit as Ax goes to 
zero. Since both A and p are constant, the continuity equation becomes 


J_ 3 (pu) . 3p _ n 

P 3x 3t 


(II-4) 
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Equation II-4 is used to describe one-dimensional mass flow in a porous 
material and is one of the three conservation equations used in 
Section D. Except for the porosity term P, Eq. II-4 is the usual con- 
tinuity equation for fluids. In all future derivations u and P are 
combined as u/P because a change in porosity only modifies the particle 
velocity in Eq. II-4. 
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B. Momentum Equation Derivation 

In the following momentum equation derivation a one- 
dimensional model is assumed. As discussed in Section II-A, the po- 
rosity is included with the particle velocity. 

The forces on the control volume fluid are illustrated in 
Fig. II-2. The conservation of momentum states that 


the time rate of change 
of momentum inside the 
control volume 


the sum of external 
forces in the x- + 
direction (pressure 
and drag) 


x-momentum 
inflow 
through the 
two end sur- 
faces . 


For the one-dimensional system shown in Fig. II-2, the momentum 
equation is 


A 3pu 


P at 


Ax * PA - PA 


x 


- FAAx + p -^r A| - p ~ A! , . 
x+Ax p2 'x p 2 x+Ax 


(II-5) 

where P is the total pressure and F is the drag force/unit volume. 
Dividing Eq. II-5 by AAx and taking the limit as Ax^O yields 


J. 
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(II-6) 


Rewrite Eq. II-6 as 


P 




P u + pu 
* t 



P 



puu 
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p uu 
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F 


(II-7) 


where, from Eq , II-4, the underlined terms equal zero. The x-direction 
momentum equation becomes 
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FIGURE 11-2 

STATIONARY CONTROL VOLUME USED IN THE DETERMINATION OF THE 
CONSERVATION OF MOMENTUM RELATION FOR POROUS MATERIALS. 
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Except for the P and F terms, Eq . II-8 is the momentum equation for 
inviscid flow in an open fluid. 

The drag term F in Eq . II-8 must be evaluated. Fulks 
et al. [29] give a derivation of their three-dimensional term which 
is equivalent to F. A different approach is used here. By using the 
following analysis F can be defined experimentally. 

Both the drag and the dc flow resistivity are used to 
define dissipation in a porous material. Whereas we have not defined F, 
we can measure the flow resistivity and mathematically relate the two 
variables. The dc flow resistivity in a porous material is determined 
by first forcing a gas through the material at various, known particle 
velocities. The pressure drop AP is measured at each particle velocity 
u/p. Once these measurements are obtained the dc flow resistivity is 
calculated by using [2] 


PR 


dc . Pap 
uL 


(II-9) 


where L is the sample length. The dc flow resistance is a function 

dc 


of particle velocity and is fit to the first order equation 


[12,19,21,55,88] 
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dc 


. u 

o + n~ 


(11-10) 
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where the coefficients a and n are fit to the dc flow resistivity data. 

The drag term F is defined by noting that if the flow is 
steady (dc) , then the first and second terms of Eq. II-8 equal zero 
[86] . The momentum equation becomes 

P = -F . (II-ID 

Divide Eq . 11-11 by the particle velocity inside the porous material 
u/P to determine the dc flow resistivity inside the material 


PR. P pv 

dc _ p ,x rF 

L u u 


(II-12) 


and 


F = 


(a + n 


.h\ H 
P/P 


(II-13) 


Both Hersh and Walker [33,34] and Zorumski and Parrott [86] have shown 
the flow resistivity to be independent of the frequency. Thus, the 
coefficients o and n are independent of any time derivatives. 
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C. Energy Equation Derivation 


In this section the internal energy equation for the fluid in 
the porous material is derived for use as the third conservation 
equation in the perturbation analysis. As discussed in Section II-A, 
the porosity is included with the particle velocity. 

The energy transfer in the control volume fluid is illus- 
trated in Fig. II-3. The conservation of energy states that 


the time rate of 
change of energy 
inside the con- 
trol volume 


net flux of 
energy across 
the control + 
volume sur- 
faces 


work/unit time heat 

done by sur- added 

face forces on + to the 
the fluid fluid . 


For the one-dimensional system shown in Fig. II-3 the energy 
equation is 


AAx 
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(11-14) 
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where e is the internal energy/unit mass and q is the heat flux/unit 
area. Dividing Eq. 11-14 by AAx and taking the limit as Ax->-0, yields 


(e+„ 2 /2P 2 ) , . [p ( e + u 2 /2P 2 ) u] _ 

, t P ’ P 


(11-15) 


Equation 11-15 is rewritten as 
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FIGURE 11-3 

STATIONARY CONTROL VOLUME USED IN THE DETERMINATION OF THE 
CONSERVATION OF ENERGY RELATION FOR POROUS MATERIAL. 
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where the total derivative is 


_D_ 3 u _3_ 

Dt 3 t P 3x 


(H-17) 


The terms in the square brackets of Eq. 11-16 equal Eq . II-4, the 
continuity equation. The terms in the curly brackets of Eq. 11-16 
equal-F in Eq. II-8, the momentum equation. The energy equation 
becomes 



q 


x,x 


+ 



(11-18) 


The fluid in the porous material is assumed to be an ideal gas, 
which has the properties 


P = RpT 


and 


(IT-19) 


de « C dT 
v 


( 11 - 20 ) 


where R-(y-l)C , y=C /C , T is the temperature, and C and C are 
v p v p v 

the specific heats. By using Eqs. 11-19 and 11-20, Eq. 11-18 is 
rewritten as 


Dt 



+ 


(y-l)Pu 

pP~ 


jc = (y-1) 

p 



( 11 - 21 ) 


We may rewrite the continuity equation (Eq. II-4) as 


l£e. + 

p Dt 


u 


j2L 


P 


= 0 


(11-22) 
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Equation 11-22 is used in Eq. 11-21 to eliminate the density variable 
such that 

DP ^ Pu 


Dt 


■f 2 • 


(11-23) 


This form of the energy equation is simplified by noting that the 

heat flux q is negligible. There are three reasons for ignoring the 

heat flux term q in Eq . 11-23. First, in open air q is usually ignored 
x x 

between 20 Hz and 20 kHz, because the relatively long wavelengths and 
short oscillation times are not conducive to heat transfer. Second, in 
a porous material the wavelength is much longer than the fiber spacing 
and the heat cannot flow along the direction of propagation without en- 
countering many fibers in a wavelength. The analysis presented in 
Appendix A shows that heat transfer effects are localized around the 
fibers of the porous material, thus further reducing q^ ^ . Finally, in 
a porous material it can also be assumed that Fu/P>>q^ ^ and Eq . 11-23 
is rewritten as 


yPu 


dp + : — i* = (v-i)U F 

Dt p ' p 


(11-24) 


If the heat transfer rate from the fibers to the air is large enough, 
the compressions of the air are almost isothermal, y~l, and 

pu 

+ -p^ = 0 . (11-25) 

Equation 11-25 is the internal energy equation used in the 
perturbation analysis that is presented in the next section. For 
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high intensity sound the rate of heat transfer increases because of 

convection and mixing in the air (see Appendix A) and we assume Y-l 

for calculation of the acoustic terms of the perturbation analysis. 

Equation 11-25 leads to the relation b =/P /p . 

o o o 
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D. Perturbation Analysis 

In the last three sections the conservation equations were de- 
rived, In this section perturbation analysis is used to solve dimension- 
less forms of the conservation equations for a single boundary condition. 
The boundary condition is a second harmonic distorted sine wave. Both 
components are of arbitrary phase and amplitude. The boundary condition 
is used in an attempt to approximate the measured signals which are des- 
cribed in Chapter IV. 

1. Perturbation of Conservation Equations 

In terms of the total derivative, the three conservation 


equations derived in Sections A, B, and C are 


1 Dp = 

P Dt p 


= 0 , 


(11-26) 


£2h + P 

pDt ,x 


' °p " n (p) 


(11-27) 


and 

YPu 

5F + -p"* = (Y-l) f F . (11-28) 

We remind the reader here that, in a porous material, we will assume 
isothermal conditions (y = 1) . For generality and because y = 1 is not 
exactly true, we use y^l to compute the perturbed wave equations 
and, in solution of these equations, set y = 1. 

The measurement results, presented in Chapter IV, show that 
the initial waveform is not a pure sinusoid. The microphone is located 
some distance from the source (See Chapter III) and the intense sound 
distorts as it travels to the microphone. At high intensities, if the 
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level of the second harmonic component is large, it can affect the 
propagation of the fundamental component [81]. Excess attenuation of 
the fundamental component may occur [16,80], but this is a third-order 
effect for an initial sinusoid and can be a second order effect for a 
second harmonic distorted sinusoid. Because of the mathematical com- 
plexities an unsuccessful attempt was made to derive a third-order solu- 
tion of an initially sinusoidal wave. 

The second harmonic distorted sinusoid boundary condition is 
used in an attempt to describe the propagation of intense sound in a 
porous material. The boundary condition is used to introduce both the 
first and second harmonics at the boundary 

p(0,t) = p T (sinu)t + bsin2wt) (11-29) 


and, since we only consider outward traveling waves, p(°°,t) = 0, where p 1 
is the fundamental component amplitude and b is the harmonic component 
relative amplitude. Because of more higher harmonics present at the 
boundary, Eq. 11-29 does not accurately describe the conditions at the 
first microphone. This inaccurate description appears to be a source 
of problems in the Chapter IV data-theory comparisons. 

The first step in the analysis is to make the equation 
dimensionless. This is done by using the following substitutions: 


* P TT u P 

A= r — , u = — , n = 

p c ’ 2 

o pc 


u 

xto , 0 

X = — , t = a)t , and e = — 
c c 


(11-30) 
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where we defined the sound speed c as neither adiabatic c nor isother- 

o 

mal b . In the solution of the equations we let c -* b . The dimension- 
o o 

less equations are 


UA AU 

A + — ^ ^ = 0 

»T P P 


5 


AU 

a 

p 


UAU 



P 


+ s 



(11-31) 


(11-32) 


and 


n 


»T 



(Y-i) 2 F 


1 


(11-33) 


where a, = a/p to, tk = r,c /p cj, and F, = F/p w. The dimensionless 
x o 1 o i o 

boundary conditions are 


H (0, t) ■ e (sin t + b sin 2t) ■ eIm(e^ T + be^ 2x ) 


(H-34) 


and II (°° , x ) = 0 . 

The solution to these nonlinear equations for a wave traveling 
in the material can be found through perturbation. Each of the three 
dependent variables are expanded in a power series in e as follows [15] 


and 


n = n 

+ 

ell. 

+ 

e 2 n 

+ 

e 3 n + 

0 


1 


2 


3 

u = u 

+ 

elJ, 

+ 

C 2 V 

+ 

e 3 U„ + 

0 


1 


2 


3 

A « A 

+ 

eA. 

+ 

e 2 A„ 

+ 

3 

e A, + 

0 


1 


2 


3 


(H-35) 
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2 

where II q - P Q /c P o = l/y, 11^ = 0, and = 1. These expansions are sub- 
stituted into both the conservation equations and the boundary condi- 
tions. The terms multiplied by like powers of c n are then combined to 
give the n^* 1 order equations. 

The first-order conservation equations are 



0 


(11-36) 


and 


4^ + n 


° 1 U 1 

, + -=-= = 0 
1 .x P 




0 . 


(11-37) 


(11-38) 


These equations are the small signal conservation equations for a rigid 
porous material. Combining Eqs. 11-37 and 11-38, to eliminate 
yields the wave equation 


(11-39) 


which has been solved by others [23,30,32,33,45,68] and is discussed in 
Chapter I. The boundary conditions are 


VO,.) = sinx 4- b sin2t and = 0 


The second-order conservation equations are 


U 2 T a i U ? 

^ + n 2 v + — 


- A ! Ul » T _ u ! Ul .y _ n i u i 

p p 2 p 2 


and 


(II-4Q) 


(11-41) 


(11-42) 
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2,T p p p p 2 


(II -A3) 


Combination of Eqs. 11-42 and 11-43 to eliminate the U 2 terms yields the 
second- order wave equation 


- n„ + a,n_ _ = - -n (u,n. ..) „ + a. u n 


2, XT 2, XX 1 2,t 


1 I.x'.t 1 1 1,X 


+Y(U l,x”l ) ,x +Y0 1 D 1 ,x"i + (n i D i,t ) ,x] 

- ji [fc/i.T + + 'Vu'.x + n i u Lx] • 


(11-44) 


The second-order boundary conditions at x = 0 and °° equal zero. The 
first order solutions are used in the second order equations and then 
the second order solutions are found. 

2. First-Order Solution 

In small signal, first-order, wave propagation, each harmonic 
component travels independently, i.e., superposition holds. As in 
Chapter I we can assume a pure sinusoid 

j (x-r.x) 

n i (x,x) = Be > (11-45) 


where T. = (B , - jot.)b /w. Substitution of Eq. 11-45 into Eq. 11-39 yields 
1 1 1 o 

the dispersion relation 


I i"V 1 - Ja i . 


or, in dimensional form, 

2 2 
? GO) 

(b 1 - ja r - ^ + j-V 

b _ b 
o o 


= 0 


(11-46) 


(11-47) 
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where we have set c=b^. Equation 11-47 can be written in terms of 
the attenuation 


a, =" 7 ^ — >/- 1 + \/l + a7 
1 v/2b c 1 

and the phase 


e =-~ 
1 VTb 


v4 + \/l +~a 


As noted in Chapter I, the phase speed c^p^ = is 

\f. 2b 


U> 


'1PH 6 


1 


Ji+ 7i 


+ o. 


(11-48) 


(11-49) 


(11-50) 


which we use to describe how fast a wave travels at each frequency. The 
phase speed and attenuation are used in Chapter IV to compare the low- 
intensity theory and experiment. 

The first-order radiation boundary condition (Eq . 11-40) is 

n i (0,t) = Im(e jT ) + bIm(e J(2T+<t,) ) . (11-51) 

The first order solution of Eq. 11-39 is 

j(t-r x ) j[2(-r-r x>++] 

n (x,t) = e +be z , (11-52) 


where = A - jo^/2 . 

The first order particle velocity function is determined by 
using Eq. 11-52 in Eq. 11-38 to obtain 


U 


1 


P 


j (x-r.jX) 

e 


bP 

r 


[2(t-r 2 x) + 4>] 


(H-53) 
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We stated in Section C that when y = l,the acoustic pressure 

2 

and density are related by a simple constant b Q . In the present case we 
find that 

Hj-A . (11-54) 


3. Second-Order Solution 

Equations 11-52 through 11-54 are used in Eq. 11-44 (set y = 1) 

to obtain 


b( 2 r 2 -r*)[ 2 n 1 +jrj(r*-r 2 )] j[x-(2r 2 -r^)x-H>] 


n 2,TT n 2 ,xx + a i n 2,t 


* 

2r r 

1 2 


(-rij - j4r 2 r 1 ) j 2 (r-r lX ) 


b(2r 2 +r 1 )(n 1 +j3r 2 r 1 ) j [3x-(2r 2 +r 1 ) x +4- D 

FT e 

1 2 


2b z (n 1 +ji6r 2 rJ) j [4(x-r 2 x)+24>3 


( r i -r i> ( -ri i +r i a f) 3( r r r x )x 


* 

2r r 

l l 


b2(r 2' r 2 )(_T1 l +r 2 CT l ) j 2(r 2" r 2 )x 


* 

2r r 

2 2 


(II-55) 


The last two terms of Eq. 11-55 are streaming terms and are ignored in 
the calculation of an acoustic solution. 
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The solution of Eq, 11-55 is obtained by first writing the 
left hand side differential operator as 


, _ 3 2 d 2 9 

L " a 2 2 1 9 t 

9T 9X 


(11-56) 


and then, by using the theory of differential operators [65], the par- 
ticular solution is found to be 


n 2 p - Qje 


j[T-(2r 2 -r 1 )x-N>] j2(x-r lX ) j[3 T -(2r 2 +r 1 ) x +<t>] 


+ Q 2 e 


+ Q 3 e 


+ v 


j[4(T-r 2 x)+2<fr] 


(H-57) 


where 


and 




b(2r 2 -r*)[2n 1 +jr^(r*-r 2 )] 

8r*r 2 (i-r*r 2 ) 


^2 = 


- AF 2 r i 
2 °i r i 


q 3 = - 


b(2r 2 +r 1 )(n 1 +j3r 2 r 1 ) 

4r i r 2 (1 - r i r 2 ) 


^4 = 


2b 2 (jn r i6r 2 r 2 ) 


(11-58) 
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Since the sum of the particular and homogeneous solutions must satisfy 
the boundary condition 


II 2 (0,t) = 0 , 


(H-59) 


then the homogeneous equation has the form 


± / — v v4 .n j(2T-2r x) 

J(T-r x +*) +R 2 + 

n 2H ' R l e 2 3 


j[3(t-r 3 x)-^] 


j[4(x-r x)+2<t’] 

+ R.e , (11-60) 

4 

where the terms are to be determined by using the boundary condition 
Eq. 11-59, The propagation parameters are Y ^ = /l - ja^ /3 and 

- /I - j a -^ / 4 . The form of the propagation parameters is determined 
by substituting the individual exponentials into the operator L of 
Eq, 11-56 and solving the dispersion relations. 

We add the homogeneous and particular solutions, Eqs. 11-60 
and 11-57, and substitute the result into the boundary condition. We 
find that -Rj=Q^ and and obtain 


/ -j( 2 r -r*)x -jr x 

n 2 (x,T) = Q 1 (^e -e 

)e^ + *Ko 2 {: 12TiX .: >2T2X )j 2 ' 

/ -j(2r +r ) x -j 3r 

+ Q 3 (e -e 

X ) e 3(3T**) +q4 ( c - J4r 2 1 < _ e - J4r 4 X ) e j(4T + 2„) 


(11-61) 
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To check Eq, 11-61 against a similar equation in lossless, 
isothermal air we used the substitution P+1-6 and took the limit as 
6 + 0. The equations were equal. 

Although Eq. 11-61 is only of second order, the equation 
yields information on four harmonic components. Two of the components 
start at nonzero amplitudes, whereas the other two start at zero ampli- 
tude. As we shall see in Chapter IV, this representation of the actual 
measurement situation is not correct, but, in order to obtain a more 
realistic representation of experimental results the boundary condition 
would also have to include the third and higher harmonic components. 
Although this latter problem is tractable, the mathematics are prohibi- 
tive for the relative information that the solution would yield. The 
solution has not been attempted. 

Figure II-4 shows a case of propagation for assumed values of 
the material properties in Eq. 11-61. Computer program PERT4PD was 
used to generate this plot and is listed in Appendix B. We note several 
phenomena illustrated on this plot. We see that the intensities of the 
higher harmonics grow with distance and then decay . The fourth harmonic 
is not shown because its level is less than 112 dB. The fundamental 
propagates in a linear fashion. The fundamental propagates in a nonlin- 
ear fashion at higher sound levels, for both e and b, and for higher 
relative nonlinearities n /a. The nonlinear propagation is exhibited by 
extra attenuation of the fundamental and more energy being transferred 
to the harmonic components. When the initial amplitudes are high 
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FIGURE 11-4 

SOUND LEVEL VERSUS DISTANCE FOR THE FIRST FOUR HARMONIC 
COMPONENTS OF A SECOND HARMONIC DISTORTED SINUSOIDAL 
WAVE PROPAGATING IN A POROUS MATERIAL. 




enough, the second harmonic level exceeds the fundamental level and the 
model fails. More discussion is presented in Chapter IV, where we find 
the measured data is not well predicted by the theory. 
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E. An Approximation to the Impedance and Absorption Properties of 
Semi-Infinite Nonlinear Porous Materials 

As discussed in Chapter I, several investigators [23,32,33, 
45,68] have determined relations for the impedance and absorption 
properties of linear porous materials. Experimentally, the impedance 
has been found to change with intensity [47,55,64,86]. In this section 
the dimensionless impedance and absorption equations are found for 
linearly behaving materials. Because acoustic impedance of porous 
materials is frequency dependent, the impedance concept is not 
easily applicable to nonlinear systems. Thus, only an estimation 
for nonlinear behavior is made for the absorption and impedance. 

The acoustic impedance of porous materials is found in 
the same manner as in Chapter I. For linear behavior the first order 
solutions U ] _ and IT^ , expressed in Eqs. 11-53 and 11-51 (b = 0) , are 
used in calculating the dimensionless impedance 


Z 

1 



= /l - 


The resistive (real) part is 


Re(Z 1 ) = 



and the reactive (imaginary) part is 


(11-62) 


(11-63) 


Im(Z 1 ) = - 



(11-64) 
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These functions are the dimensionless counterparts of Eqs. 1-8 and 
1-9. 

Once the propagation properties of a high-intensity sound 
in a fibrous porous material are known, it is advantageous to determine 
the impedance that an intense incident wave encounters at the material 
surface. The acoustic resistance that an intense sound encounters is 
influenced by the material nonlinearity and the amplitude of the sound. 
If an attempt is made to solve for the ratios II /U or the 

results are complicated functions of the different harmonics and their 
relative amplitudes. Although a solution might be obtained by using 
the perturbation solutions, it is expected that the result would be 
far more complicated than the perturbation analysis and not very useful. 
A simpler approximation is presented here. 

The following approximation is good only at the air/surface- 
interface of a semi- inf inite nonlinear porous material. For a sound 
traveling in the material the perturbation solutions or one of the 
approximations in Section II-F should be used. To avoid the problem of 
reflection from the back of the sample, we assume that aL >> 1. 

We employ a simple substitution in the approximation. In 
place of the linear dc flow resistivity a^, we use the first order 
relationship for the nonlinear dc flow resisitivty o ^ + . This 

substitution is similar to that used by Ingard [37] for high-intensity 
impedance of Helmholtz resonators. We use the quantity e//2 as an 
acoustic replacement quantity of the dimensionless dc particle 
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velocity u^ c /b Q . We do not want any time or distance dependency 
and work only with the appropriate magnitudes. The substitution 
of this relationship into Eq . 11-62 yields the dimensionless nonlinear 
impedance 

Z N = \/l - j (<?! + en 1 /P ,/, 2) » (11-65) 

with the resistive and reactive parts becoming 


Re(Z 


) ■ Vl + 7l + (o +en,/p/ 2 ) 2 

sjl 1 1 


and 


(11-66) 


Im(Z ) 
N 


_=L 

& 


y-i+Ti + (a^enj^/p/2 )' 


(11-67) 


For a specific material, the above functions may be plotted 
with respect to the parameters of frequency and acoustic Mach number of 
the incident sound at the surface. In Fig. II-5 the normal specific 
acoustic impedance encountered by a 1 kHz wave at the surface of a 
semi-infinite, fibrous porous material is plotted with respect to 
the sound pressure level. The magnitudes of both the resistive and 
reactive parts of the impedance increase with amplitude. These re- 
sults Indicate that both the amplitude and phase change of the 
reflected wave vary as the amplitude increases. In the analysis 
review presented in Section I-C-2 we found that the reflection process 
appears to follow the small signal laws up to about 170 dB. In the 
test procedure described in Section III-E-2 we find that the results 
are reliable up to a sound level of about 165 dB. The comparison 
between the theory and data should be appropriate. 
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FIGURE II 5 

ACOUSTIC IMPEDANCE VERSUS SOUND PRESSURE LEVEL FOR A 
SEMI-INFINITE POROUS MATERIAL. 



The effect of the increased resistance and reactance on the 
ability of the sound to enter the porous material is examined by 
studying the absorption coefficient. We can rewrite Eq. 1-11 in the 
nonlinear dimensionless form 


a N = 1 ' 


Z N " 1 


Z + 1 

N 


( 11 - 68 ) 


The nonlinear absorption coefficient at 1 kHz versus sound 
pressure level is plotted in Fig. 11—6 for the same material properties 
as used in Fig. II-5. The absorption coefficient decreases with 
increases in the sound level. When the resistance increases, the 
absorption coefficient is expected to decrease. Since the flow resis- 
tance increases with sound level, it is harder for the high-intensity 
sound to enter the material than for the same frequency low- intensity 
sounds. 

The increased impedance has two consequences in the absorp- 
tion of sound. First, the high- intensity sound cannot enter the 
absorbing material as easily as low- intensity sound. Second, once 
the high- intensity sound does enter the material it is more rapidly 
dissipated than the low- intensity sound. When the level is reduced to 
low intensities, the sound asymptotically approaches linear propaga- 
tion. The latter effect is noted in the experimental results (Section 
IV-B-2) and in the Section II-F approximations. 
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FIGURE 11-6 

ABSORPTION COEFFICIENT VERSUS SOUND PRESSURE LEVEL 
FOR A SEMI-INFINITE POROUS MATERIAL. 
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F. Approximations to the Attenuation of an Intense Sinusoid 
Propagating in Nonlinear Porous Materials 

In this section we present two models for use in approximating 
the attenuation of intense sinusoids in porous materials. The linear 
impedance and propagation parameter V ^ are intimately related by 
Eq. 11-62. Because of this relation, we can use the approximation of 
Eq. 11-65 to model how the propagation parameter of an initially 
sinuoisdal wave changes with intensity. We may also analyze the 
propagation of an initially sinusoidal wave by modeling how energy is 
lost from the fundamental component as a function of sound level. The 
mathematics of each model is simple and the models are easy to use. The 
impedance model is presented first. 

1. Impedance Model 

We can write the dimensionless propagation parameter as 

r i = Bi - jA x , (11-69) 

where the dimensionless attenuation A^ = a^b^/tu is given by Eq. 11-63 
and the dimensionless wave number P is given by Eq. 11-64. 

The resulting equations are analogous to Eqs. 1-4 and 1-5. 

Again, the substitution of a^+cri^/P/2 for is used to 
replace the linear equations with 


Im(r N ) = ^ = - — yfl + x/l + ( a 1 + e n 1 / P >^') 2 (11-70) 

/l 


for the attenuation and 


Re (r 1 ) = b n 


_1_ 

Jl 


\J l + Jl+(o 1 + er il /P/2)^ 


(H-71) 
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for the wave number. 


The dimensionless phase speed C j^/b Q = 1/B^ 


is 


^phN ^ 2 

b ■ J 7 = «!- 72 > 

° y 1 + vl + (c^ + zt\^IV7i ) 

We conclude that both the phase speed and attenuation change 

with intensity. Because of the interaction of the harmonic components, 

measurement of C , % ,/h is difficult; one cannot differentiate between 
phN o 

actual phase speed reduction and harmonic interaction. Of course, 
the effects of the interaction may be to slow as well as attenuate the 
fundamental component. In some of the experimental results we measured 
an apparent slowing of the fundamental, but, because of the limited 
data, we do not present the results in this study. 

In Fig. II-7 the amplitudes of waves propagating in porous 
materials are plotted versus distance. The solid line indicates the 
effect of linear attenuation. The dashed line indicates the effect 
of nonlinear attenuation (Eq. 11-70). We have used exp[/A^dy] to 
determine the amplitude attenuation. The attenuation is more rapid 
where the intensity is highest. The difference in the attenuation is 
called the excess attenuation. In Chapter IV the approximation, 

Eq. 11-70, is compared to the perturbation solutions and data. 

2. Amplitude Attenuation and Saturation Model 

Excess attenuation of the fundamental component of an initial- 
ly sinusoidal wave may also be modeled by a simple rate equation, which 
includes the decay of the fundamental component due to nonlinear effects 
as well as the direct decay due to small signal dissipation. We base 
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FIGURE II 7 

SOUND PRESSURE LEVEL VERSUS DISTANCE FOR AN INTENSE 
SINUSOID PROPAGATING IN A POROUS MATERIAL. 




our model on the concept that, at high intensities, the fundamental 
loses energy in a consistent manner, regardless of the medium. 

In Chapter I we found that the rate of harmonic component 
generation depends on the amplitude of the original sinusoid and on the 
thermoviscous attenuation of the fluid medium [l6]. The harmonics are 
generated at the expense of the fundamental. If the rate of harmonic 
generation is high enough and the attenuation low enough, a shock will 
form. The presence of the shocks Increases dissipation because losses 
are large at the shocks. Energy Is lost by all the harmonics. Thus, 
we find that the greater the amplitude of the original sinusoid, the 
quicker the fundamental loses energy. 

If we measure the sound level of an intense fundamental 
component at some remote distance from the source, we will find a 
nonlinear relation between the source level and remote level. Because 
the overall decay of the fundamental depends on the source level, the 
remote sound level depends on the source level, the small signal attenu- 
ation, the nonlinearity of the medium, and the distance from the source. 
A* limit in the remote sound level, which is called the saturation level, 
is a function of the above variables. The excess attenuation, linear 
region, approach to saturation, and saturation level are all illustrated 
in the sketch of the amplitude response plot in Fig. II-8, In Chapter 
IV we present data showing the approach to saturation. We model the 
saturation effects in the following analysis. 

The excess attenuation that the fundamental component of a 
sinusoidal wave encounters In a porous material may be modeled by 
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FIGURE 11-8 

AMPLITUDE RESPONSE CURVE FOR AN INTENSE SINUSOID 
PROPAGATING IN A POROUS MATERIAL. 
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considering the two factors causing attenuation of the fundamental 
component. Others [4,57,76,80,81] have used this method to explain 
saturation effects in plane waves traveling in an air-filled tube. 
Ordinary small-signal decay dominates attenuation at low intensities 
and is described by 


Pj_ = P Q e “ X , (H-73) 

where p^ is the pressure amplitude of the fundamental component, 
is the initial pressure amplitude of the fundamental component, and 
x is the distance from the source. The decay rate equation associated 
with Eq. 11-73 is 


dp i 

d T = " ap i 


(11-74) 


Nonlinear effects dominate the attenuation at high intensities 
and loss of energy from the fundamental harmonic component depends on 
intensity. Webster and Blackstock [ 81 ] made measurements and confirmed 
that the pressure amplitude of the fundamental component in the sawtooth 
region is 


2 P. 


(H-75) 


P 1 ~ 1+ o* 

where a* > 3. Their assumed boundary condition is that the wave starts 
as a sawtooth. Since a* = x/x and a* > 3 the pressure relation is 
approximately p-^l/x and the decay rate equation associated with 
Eq. 11-75 is 


dp x 


dx 


a 


(11-76) 
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Our measurement results of intense sound propagating in a 
porous material have shown that excess attenuation exists. We also 
found that shocks do not form in a porous material. In spite of the 
lack of shock formation, the approach to saturation appears to be the 
same as that found in air. We propose that the high-intensity decay 
rate equation in a porous material is 


dp 1 


dx 



(11-77) 


where T is an unknown coefficient dependent upon the nonlinearity 
of the medium and the fundamental component frequency. The coefficient 
T may be defined from either saturation or propagation tests. In 
Chapter IV we define T in the saturation tests. 

Let us assume that the two rates of decay, Eqs. 11-74 and 
II“77, may be added to obtain the overall decay at any pressure p . 

We obtain the relationship 


dx 



(11-78) 


whose solution satisfying 


the boundary condition 


p i 


p at x=0 is 
o 


P e 

Pi = ~ — ITi • (11-79) 

l+P o (l-e X )/aT 

Equation 11-79 is used in Chapter IV to describe both saturation and 
propagation of intense sinuoidal waves. 

Equation 11-79 is generally useful. At low source amplitudes 
P G ^ aT ^q* 11-79 reduces to Eq . 11-73. At high source amplitudes 
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>> aT and Eq. 11-79 may be used to calculate the saturation level. 
As p Q increases without limit p^ reaches the saturation pressure and 


Eq. 11-79 becomes independent of p^. The saturation pressure is 


aTe 


-ax 


Is 


~ ax \ 
(1 - e ) 


(11-80) 


At low attenuations a 0 Eq . 11-79 reduces to 

From the above analysis, we conclude that the three variables 

T, a, and p help us define how a sound wave attenuates. In Chapter I 

2 

we found that the Gol f dberg number G=$ , ek/a or B'p k/ap c , is used 

o o o 

to define the relative effects of attenuation and nonliner distortion. 

We conclude that, for a porous material, a Gol'dberg type number may be 
defined from Eq. 11-79 as 

G = ^ . ( 11 - 81 ) 

We can use Eq. 11-81 to describe when nonlinear effects become important 
to the propagation of intense sounds in a porous material. As 
nonlinear effects become dominant and as G^O attenuation effects 
dominate. In Chapter IV we determine T and relate T to the relative 
nonlinearity of the material r\/o and the frequency f. 


82 



G. Summary 


In this chapter the propagation of sound In and the impedance 
of porous materials have been modeled. At low intensities the theory 
presented here agrees with work done by others. At high intensities 
harmonic components of the fundamental are generated and the attenuation 
of the fundamental component is increased. The Increase in flow 
resistance at high intensities is found to effectively increase the 
impedance of the materials. The impedance increase makes it more 
difficult for a sound to enter and propagate in the material. 

In the next chapter, Chapter III, the methods for testing 
and evaluating porous materials are discussed. In Chapter IV the 
results of Chapter II are compared to data obtained from the experi- 
ments and conclusions about the theoretical methods are discussed. 



CHAPTER III 


EXPERIMENTAL METHODS 

This chapter contains a description of the experimental methods 
and physical apparatus that were used to determine the properties of the 
bulk porous materials. The chapter is divided into the following 
sections : 

A. Porosity 

B. Material Structures 

C. DC Plow Resistivity 

D. Acoustic Propagation Parameters 

E. Specific Normal Acoustic Impedance 

F. Summary 
A. Porosity 

The porosities of the various materials were determined by 
using the following procedure. The density of each sample was deter- 
mined by carefully measuring the sample and weighing it on an analytic 
balance. The density of the skeleton material was determined from the 
literature or the manufacturer. The porosity was determined from 



where m and V are the measured mass and volume, respectively, of the 
sample and D is the density of the solid material. The porosity ranges 
of the various materials used in the present study and some densities 
for the solids of acoustical materials are listed in Table III-l. 
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TABLE III-l 


A LIST OF THE POROSITY RANGES FOR THE SAMPLES USED IN THE PRESENT 
STUDY AND DENSITIES OF VARIOUS ACOUSTICAL MATERIALS 


Material 

Density Range 
D,kg/m 3 

Porosity Range 

P 

Kevlar 29 aramid [ 90] 

(an aromatic polyamide) 

1439 ± 30 

0.985-0.809 

Polyimides [ 5 o] 

1439 

— 

Scottfelt 900-Z [89] 

1153 ± 32 

0.942-0.853 

Polyurethane Plastic [ 5 0] 

1107 - 1246 

0.971 

Glass (minimal binder) [ 5 0] 

2518 - 2601 

0.982 

Aluminum [ 5 0] 

2600 - 2900 

— 

Steel [50] 

7400 - 7800 

— 
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Beranek [6] devised an instrument to measure the porosity di- 
rectly. A device of this type was built for the present study. The 
device is not as accurate and the measurements are not as repeatable as 
the weight method described above. 
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B. Material Structures 

The material structures were determined by photographing the 
materials with a camera-bellows combination. Fiber size and material 
structure are useful in determining the material properties. Photos 
of the Kevlar 29 (manufactured by du Pont de Nemours & Co., Wilmington, 
DE) and the Johns-Manville 1000 (J-M 1000) fiberglass (supplied by 
Johns-Manville) are shown in Fig. III-l, of the Scottfelt 900-Z-2 
(supplied by Scott Paper Co., Foam Division, Chester, PA) and the 
Blachford Acoustical Foam (BAF) (supplied by H. L. Blachford, Inc., 
Corona, CA) in Fig. III-2. 

We see that the Kevlar 29 fibers are larger (12 pm) than the 
glass fibers (10 pm). The structural members of the two foams are 
larger than the Kevlar fibers. The Scottfelt has an effective diameter 
£ 40 pm (±15 pm) and the BAF effective diameter ranges between 30 and 
140 pm (we use 75 pm in later calculations). The BAF has "window 
panes," which partially close off some of the cells of the foam. The 
BAF is a partially reticulated foam, whereas the Scottfelt is a fully 
reticulated foam, i.e., the Scottfelt has no "window panes." 
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A. SCOTTFELT 900-Z-2 



B. 


BLACHFORD ACOUSTICAL FOAM 




I mm 

1 


FIGURE III-2 

PHOTOGRAPHS OF SCOTTFELT 900-Z-2 AND 
BLACHFORD ACOUSTICAL FOAM 
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C. DC Flow Resistivity 

In Chapters I and II the dc flow resistivity a of a porous ma- 
terial was shown to be important in determining the acoustical properties 
of a porous bulk absorber. In Chapter II, Eqs. II-9 and 11-10, we found 
that the flow resistivity can be measured in dc flow tests. The 
measurement method is straightforward and follows the ASTM standard [2]. 
The details of the particular measurement system used here are pre- 
sented, 

A block diagram of the measurement system is shown in 
Fig. III-3. The porous samples were carefully cut and placed in the tube, 
which was made of clear Plexiglas. The clear plastic enables the ex- 
perimenter to accurately measure the material length at any flow rate. 
The in place length measurement is important because some materials 
tend to compress with an increase in flow rate. If the compression 
goes unnoticed, an anomolous measure of the flow resistivity is ob- 
tained. Most samples were 7.7 cm in length. The sample was held in 
place between a fixed screen and a moveable screen. Each screen was 
made of standard h in. hardware cloth. The sample was located so that 
its ends were at least 5 cm from the intake and exhaust ports in the 
tube end pieces. There was then a fairly even flow across the surface 
of the sample. Great care was taken to avoid flanking, that is, flow 
around the outside of the sample. If flanking was noticed, it was 
stopped by inserting thin plastic sheets between the tube wall and the 
material. Corrections for the presence of the plastic sheets were 
then made in calculating the flow resistivity. 
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BLOCK DIAGRAM OF THE DC FLOW RESISTIVITY 
MEASUREMENT SYSTEM. 

j air flow meters volume flow rates (scf/hr) 

DWYER RMA-1 
DWYER RMA-3 
DWYER RMB-50 
GILMONT B31 2 
DWYER RMB-55 

| MANOMETERS 

i DWYER 1 INCH GAUGE 

WEHLING 3 INCH GAUGE 

| 

j UTUBE 15 INCH 

I 

{ r 

| 

FIGURE 111-3 

BLOCK DIAGRAM OF THE DC FLOW RESISTIVITY MEASUREMENT 
SYSTEM AND A LIST OF DEVICES USED. 
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1.0 - 10 
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40 - 400 
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Five flow meters were used to measure the velocity of the 

flow. The meter types and ranges are listed on Fig. III-3. The meter 

ranges overlapped and the total range of particle velocity measurement 
-4 

was from 7*10 m/s to 1.6 m/s. Flow rate measurements at the lower 
end of each flowmeter range tended to be in error by up to ±12%. In 
the least squares fit of the data this error was minimized because of 
the many measurements taken and because of the data overlap. 

Three manometers were used to measure the pressure drop 
across the sample. The manometer types and ranges are listed in 
Fig. III-3. The range of pressure measurement was from 5 Pa (5’10~ 5 atm) 
to 7700 Pa (0.076 atm). 
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D. Propagation Parameters 

First, attenuation a (nepers/meter) and phase 6 (radians/ 
meter) were measured in a square duct at low sound intensity for the 
various porous materials used in this study. From the phase measure- 
ments the phase speed was calculated. A simple two microphone method 
was used. Several different types of low-intensity sound signals were 
used and gave comparable results. The system and measurements are dis- 
cussed in subsection 1 below. The high amplitude sound testing was 
done with the same experimental system, but only pulsed sinusoids were 
used as signals. These measurements are discussed in subsection 2 be- 
low. 

1. Propagation at Low Intensities 

Measurements of the small-signal propagation parameters were 
made with the material to be tested in a traveling wave tube, A block 
diagram of the propagation parameter measurement system is shown in 
Fig. III-4. See Appendix B for detailed drawings of the system. The 
Brttel and Kjaer type 4136 microphones (k in. diameter) were placed in 
holes in the tube wall. The microphones were recessed (~ 4 mm) from 
the material surface. The remote microphone was vibration isolated 
(see Appendix B) . Vibration isolation for the source microphone was 
found to be unnecessary. Measurements were taken at several microphone 
separations d (see Fig. III-4). Unused microphone holes were plugged 
and sealed. The material length was such that reflections from the far 
end of the material were unimportant [a(L-d)>> l]. The passband of the 
James B. Lansing 375-H (aluminum diaphragm) acoustic driver was 
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WIRE SCREEN LOCATED 9 cm FROM THE DIAPHRAGM. 


FIGURE 111-4 

BLOCK DIAGRAM OF THE TRAVELING WAVE TUBE SYSTEM USED 
IN THE PROPAGATION PARAMETER AND SATURATION TESTS. 
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approximately 0.25-8 kHz. The upper limit of the useful frequency range 
of each tube is the first cross mode cut-on frequency in the material. 
The cut-on frequency depends on the material and is between 5 and 9 kHz. 

The propagation parameters a and 8 were measured by using 
several excitation signals, broadband noise, narrowband noise, sine 
sweep, pure tone, and tone bursts (the tone bursts are described in the 
next subsection) . In each case the fast Fourier transform (FFT) spec- 
trum analyzer was used to display the waveforms, spectra, and transfer 
function. The different excitation signals gave equivalent results and 
the transfer function of the broadband noise signal between the micro- 
phones was used to determine most of the small signal propagation 
parameters . 

2. Propagation and Saturation at High Intensities 

The experimental setup in the high-amplitude sound propagation 
testing was the same as described in Section B-l. Because of the 
limited steady-state power handling capabilities of the driver (= 25 W) , 
the high-amplitude sounds could only be achieved by using unfiltered 
pulsed sinusoids. The maximum pressure amplitude was limited 
by the power output of the amplifier (200 W at 70 Vrms) . It was 
felt that the drivers would be capable of handling this power as long 
as the duty cycle was very short [64], It turned out that one driver 
failed after several thousand pulses and a second started to generate 
unwanted harmonics after several thousand pulses. In each case the 
voice coil separated from the diaphragm. In the first case the coil 
distorted, rubbed the magnet, and shorted out. The sound pressure 
levels of the fundamental component obtained were up to 173 dB at 1 kHz, 



I 


170 dB at 2 kHz, and 172 dB at 3 kHz. These values were obtained by 
driving the amplifier up to its clipping limit (100 V peak ampli-) 
tude) . The pulse length was 0.004 to 0.005 seconds, and at least five 
pulses were averaged by the Fourier transform analyzer. The second 
harmonic distortion of the microphone system is s2.5% at 173 dB SPL, 
i.e., the second harmonic component introduced by the microphone system 
is 141 dB [20]. Above 170 dB the measured second harmonic distortion of 
the waveforms ranged from 7-25% at the source microphone. 

Two types of tests were made. The change of the amplitude 
of the harmonic components with distance was measured. During the 
same tests the distance was fixed and the change of the amplitude of 
the harmonic components with increase in the source level was measured. 
At each frequency the sound levels were started at 120 dB and increased 
in 5 or 10 dB steps until the maximum level was obtained. The latter 
study shows the approach to saturation of a wave traveling in the 
material. The phase change with amplitude between the components was 
also measured. The test results indicate how sounds propagate in 
porous materials with respect to distance, frequency, amplitude, and 
material attenuation and nonlinearity. An interesting point is made 
here about the waveforms: Even at very high intensities very little 

distortion of the waveform was seen; most waveforms were sinusoids or 
slightly distorted sinusoids. This slight distortion may be inferred 
from the spectrum measurements presented in Chapter IV. 
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E. Specific Normal Acoustic Impedance ? 

Specific normal acoustic impedance was measured by using a 
standard standing wave impedance tube for the low-intensity tests [l] 
and single cycle sine wave tone bursts for the high-intensity tests [47]. 

The low-intensity tests are discussed first. 

1. Impedance at Low Intensities ‘ 

Specific normal acoustic impedance at low sound amplitudes 
was measured with two standing wave tube systems, one for low frequen- 
cies and one for high frequencies. The block diagram of the basic 

measurement system is shown in Fig. III-5. The low-frequency (0.1- f 

1.3 kHz) tube had been constructed years earlier and is described in 
Ref. 71. The high-frequency (0. 8-3.9 kHz) tube was designed for this 
project. Construction drawings for it are included in Appendix B. 

The test method follows ANSI/ASTM C384-77[l]. The equipment 
was calibrated by using a tight-fitting aluminum plug (5 cm long) in the 
high-frequency tube and a thick steel end-cap, sealed around the edges, 
at the end of the low-frequency tube. In each case the termination was 
assumed to be of infinite impedance and the distance between the micro- 
phone and the termination was acoustically measured at each test fre- [ 

quency. The calibration results were used to correct for the environ- 
mental effects (such as tube wall attenuation, effective microphone 
position, etc.) in calculating the material impedance. In the high- 

frequency tube the probe position can be located to within ±0.1 mm and j 

i 

test results are very repeatable. In the low-frequency tube the micro- 
phone location can be determined to within ±4 mm and the error affected 
the measurement accuracy. 
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FIGURE 111-5 

BLOCK DIAGRAM OF THE STANblNG WAVE IMPEDANCE TUBE 
SYSTEM USED IN THE MEASUREMENT OF THE IMPEDANCE 
OF FINITE THICKNESS MATERIALS. 
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Materials of various types and thicknesses were measured in 
each tube, but the most common thickness was approximately 7.7 cm. (For 
checking the reliability of the single-cycle pulse measurement method, 
described next, a 30.7 cm thick sample was used in the high-frequency 
tube . ) 

2. Impedance at High Intensities 

The single cycle tone burst test used here is a simple and, 
apparently, unique method for measuring the impedance of bulk acoustic 
materials at high intensities. There are several reasons for using a 
tone burst to measure the high-intensity impedance. We have already 
discussed several reasons in Section I-C-2. 

As sound intensity increases in a standing wave tube nonlinear 
effects become evident. Distortion and shock formation cause the method 
to fail because the waveshape and amplitude change as the wave travels 
toward the sample and again after reflection. The properties of intense 
standing waves are much harder to predict than those of intense 
traveling waves. If the sample behaves nonlinearly, additional wave 
distortion is introduced at the material surface and separation of the 
nonlinear effects becomes extremely difficult. 

We can limit the noticeable nonlinear effects on the 
traveling wave by using a tone burst. As we shorten the tone burst 
less travel distance is necessary to resolve the incident and reflected 
waves at the microphone. As the amplitude increases, the shock for- 
mation distance is reduced, according to x= l/B T ek, and the waveform 
distorts at a faster rate [ 1 6 ] . When x< x, where x is the travel 
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distance, a shock is formed in the calibration signal and the 
measurement method is no longer reliable. 

The assumption on which the test method is based is that, as 
long as a shock is not formed, the attenuation of a wave as it travels 
down the tube and is reflected from the termination will proceed in an 
almost linear fashion [16]. The calibration tests at the various 
amplitudes and frequencies can be used to cancel out any small excess 
attenuation encountered by the fundamental in the production of the 
higher harmonic components. 

A block diagram of the measurement system is shown in 
Fig. III-6. The mirophone was placed to minimize the travel distance 
of a 1 kHz single cycle tone burst. The useful frequency range is 
1 kHz (set by tube length) to 4 kHz (set by cross-modes) , The sound 
level limits are dependent on tube length and frequency. When 
x=x=96.4 cm, the sound level upper limits are, for 1 kHz, 167 dB and, 
for 4 kHz, 155 dB, 

The polarity of the signal is such that a 165 dB, 1 kHz 
single cycle sine wave, shown in Fig, III-7.A, is produced. The fre- 
quency spectrum of the signal is shown in the right column. In 

Fig. III-7.B an ideal single cycle sine wave and its spectrum are shown. 
By comparing the waveforms in Figs. III-7.A and III-7.B, we find that 
the actual signal is asymmetric and contains an extra "tail.” The 
"tail" is formed by high frequency cross-modes that are generated in 
the tube and travel at different phase speeds than the tone burst. 

Both the asymmetry and "tail" affect the frequency spectrum by adding 

more low-frequency and high-frequency energy. The high-frequency lobes 
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FIGURE 1 1 1-6 

BLOCK DIAGRAM OF THE TRAVELING WAVE IMPEDANCE TUBE 
SYSTEM USED IN THE LOW- AND HIGH-INTENSITY SOUND 
IMPEDANCE TESTS. 
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IDEAL AND MEASURED WAVEFORMS AND SPECTRA OF A SINGLE 
CYCLE SINE WAVE AND OF A SINGLE CYCLE SAWTOOTH WAVE. 
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of the actual signal are not as uniform as for the spectrum of the ideal 
wave. In addition, the peak frequency is approximately 100 Hz lower 
than for the spectrum of the ideal wave. If the analysis is done in 
the frequency domain, then this extra energy can cause severe problems 
which could lead to erroneous analysis. We restrict our analysis to 
measuring the peak-to-peak amplitudes of the waves and the time delay. 

The waveform polarity was chosen such that, as the wave 

travels, a shock forms at the center. If the wave is symmetric, all 

three axis intersection points will travel at the adiabatic sound 

speed c and the frequency spectrum peak will not shift. The waveforms 
r o 

are asymmetric and after a shock forms the relative time delay is un- 
certain. If the waveform polarity was reversed, shocks would form on 
each end of the wave and the frequency spectrum would shift as a 
function of amplitude. 

In Fig. III-7.C the incident wave of Fig. III-7.A is shown 
after it was reflected off an aluminum termination. The waveform has 
almost formed a shock. The associated frequency spectrum shows that 
energy has been shifted to the higher frequencies. The wave is on its 
way to forming a single cycle sawtooth waveform. In Fig. III-7.D an 
ideal single cycle sawtooth waveform and spectrum are shown. We see 
that the spectrum in Fig. II1-7.C shows that the reflected wave is 
partway between the ideal waves shown in Figs. III-7.B and III-7.D. 

For higher frequency tone bursts the spectra are influenced 
more and more by the "tail." The resulting frequency spectra tend to 
"smear” and a localized frequency peak is not definable. If the "tail" 
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is cut off, the peak is more localized, but other high frequency com- 
ponents are introduced by the cutting operation. 

In this testing program we measured (1) the time delay be- 
tween the incident and reflected wave-center zero crossings and (2) the 
change in peak-to-peak amplitude. The ANSI/ASTM C384-77 test method 
[l] was modified to deal with traveling waves. The amplitude change was 
used to determine the material absorption coefficient and the delay time 
was used to determine the phase shift. The impedance was calculated 
from a combination of the amplitude change and the phase shift. The 
system was calibrated by replacing the material with a tight-fitting 
aluminum plug and by assuming that the measurement results indicate an 
infinite impedance. In the experiment, shock formation caused excess 
attenuation of the traveling wave. Past the shock formation distance 
calibration canceled some of the excess attenuation, but not all of it. 
In our case, asymmetry of the waveform (the positive pressure magnitude 
is greater than the negative pressure magnitude) caused the shock to 
propagate faster than the speed of sound and the associated phase shift 
could not be computed past x. This waveform asymmetry is caused by 
inertia of the driver diaphragm. 

The advantages of this measurement method over other measure- 
ment methods [59,64] are that the travel distance is short and the tone 
burst does not lengthen as it travels [ 16] • In addition, as the shock 
forms, the tone burst spectrum changes only a moderate amount. This 
spectrum change is illustrated in Fig. III-7 for both the ideal and 
measured signals [l8]. Finally, the material needs only to be long 
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enough that a wave reflected from the material termination does not 

interfere with the wave reflected from the material surface. 

The disadvantages of this measurement method are that it is 

amplitude limited and the waveform tends to be asymmetric. In addition, 

the spectrum is not that of a pure tone, but narrowband (see Fig. III-7, 

where, ideally, the -3 dB bandwidth is 0.74 f and -6 dB bandwidth 

1.1 f , where f Is the nominal frequency peak in the spectrum). 

P P 
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F. Summary 


The various measurement systems and methods for the deter- 
mination of porous material properties have been described. The topics 
covered have been on determination of porosity, material structures, 
dc flow resistivity, propagation parameters, and impedance. 

In Chapter IV the measured material properties are compared to 
theoretical results from Chapters I and II. The dc flow resistivity 
results are found to depend on the porosity and particle velocity in 
the expected manner. The material structure has definite effects on 
the acoustic measurement results. The low- intensity and high-intensity 
results are illustrative of the topics discussed in the previous chapters 
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CHAPTER IV 


EXPERIMENTAL RESULTS 

In this chapter data from experiments are compared to results 
of the theory developed in Chapters I and II. The objectives of this 
chapter are to illustrate the acoustical properties of porous materials 
with changing sound level and to determine how accurately the experi- 
mental results can be predicted by using the theoretical models. A 
variety of materials was used in the various tests and, where reason- 
able, representative results from the tests for each material are 
presented. Generalizations are made for various materials and limita- 
tions of the theories are discussed. 

In Section A we present measured and predicted results and 
discuss (1) the linear and nonlinear dc flow resistivity coefficients, 
(2) the applicability of the previously derived dependence of flow 
resistivity on porosity, and (3) the applicability of the modified 
Reynolds number in predicting nonlinearity threshold. The results of 
Section A are used in the calculations of Sections B and C. 

In Section B we present measured and predicted results and 
discuss (1) the small signal attenuation and phase speed, (2) excess 
attenuation and the approach to saturation, and (3) high- intensity 
propagation . 


107 



In Section C we present measured and predicted results and 
discuss (1) the small signal impedance for finite and semi-infinite 


porous materials, and (2) the impedance change of semi-infinite 
materials with sound intensity. 

In Section D we summarize the results of this chapter. 
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A. DC Flow Resistivity 

The physical measurement apparatus and method used to 
determine the dc flow resistivity of porous materials is described in 
Section III-A. After first presenting the flow resistivity data 
and determining the nonlinearity, we then discuss how well the empiri- 
cal models fit the data. 

1, Flow Resistivity Determination 

Flow resistivity data were taken in a particle velocity range 

-4 

of 2*10 -2 m/s. From a least squares fit of the data the coeffi- 

cients a and q were calculated; see Eq. 11-10. 

In Fig. IV- 1 the measured dc flow resistivities for five 
porosities of Kevlar 29 are plotted versus dc particle velocity. 

Each sample was 7.7 cm in length. For the three highest porosities, 
at particle velocities above 1 m/s, the data have been corrected for 
the small length compression that occurred. The compression had no 
measureable effect on the porosity, but noticeably affected the 
resistivity because the resistivity is inversely proportional to the 
material length. This material exhibits linear flow resistivity for 
most of the test velocity range. As discussed in Chapter I, all porous 
materials are expected to exhibit nonlinear behavior for particle 
velocities above some threshold value. The Kevlar 29 fiber diameter 
is 12 ym and we use Eq . 1-27 to calculate the nonlinearity threshold 

(R =1) to occur at 0.099 m/s for p= 0.980 and at 0.40 m/s for 
m 

p= 0.918. The theoretical thresholds are indicated by the vertical 
bars on the figures. In comparing the predicted threshold particle 
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FIGURE IV-1 

MEASURED DC FLOW RESISTIVITY VERSUS PARTICLE VELOCITY 
FOR FIVE POROSITIES OF BATTED KEVLAR 29. 




velocities with the points at which the curves actually start to turn 
up we see that the nonlinearity threshold is not precisely predicted by 
the modified Reynolds number at the lower porosities. The measured 
deviation from linear behavior is 2% for p= 0.980 at 0.099 m/s and 
7% for P= 0.918 at 0.40 m/s. Though slightly inconsistent, the 
predicted and measured thresholds are in general agreement. 

In Fig. IV-2 the measured dc flow resistivites for three 

different Scottfelt samples are plotted versus dc particle velocity. 

Each sample was approximately 7.7 cm in length. These materials did 

not compress with increase in particle velocity. Scottfelt exhibits 

linear flow resistivity behavior over as wide a particle velocity 

region as Kevlar. Even though the Scottfelt is not fibrous, we assume 

the effective "fiber" diameter of the foam material to be 40 ym 

(Fig. III-2). We calculate the nonlinear threshold to occur at 0.086 

m/s (3% deviation) for 900-Z-2 Scottfelt (P= 0.942) to 0.218 m/s 

(7% deviation) for 900-Z-6 Scottfelt (P= 0.852). As with the Kevlar 

the estimate based on R * 1 is reasonably good. 

m 

In Fig. IV-3 the measured flow resistivities for two pieces 
of Blachford Acoustical Foam (BAF) , Johns-Manville 1000 (J-M 1000) 
fiberglass, and Globe-Albany (G-A) needled and felted Kevlar 29 are 
plotted versus dc particle velocity. The two fibrous materials, 

Kevlar 29 and fiberglass, are linear over most of the particle velocity 
range. The BAF exhibits nonlinear behavior at a relatively low 
particle velocity. Although the two BAF samples were taken from areas 
within a half meter of each other in the same 7.6 cm thick sheet, 
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FIGURE IV-2 

MEASURED DC FLOW RESISTIVITY VERSUS PARTICLE VELOCITY 
FOR THREE POROSITIES OF SCOTTFELT. 
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FIGURE IV-3 

MEASURED DC FLOW RESISTIVITY VERSUS PARTICLE VELOCITY FOR 
GLOBE-ALBANY KEVLAR 29, JOHNS-MANVILLE 1000 FIBERGLASS, 
AND BLACHFORD ACOUSTICAL FOAM. 





the flow resistivities of the two samples were different by a factor 
of 2-4. This result illustrates that, for foams, variations 
that are not discernable can drastically influence the flow resistivity 
properties. In turn, the acoustical properties may also be expected 
to vary drastically. The porosities of each of the many BAF samples 
tested were measured to be 0.971 ±0.001. The small porosity range and 
large flow resistivity range indicate that structural factors are 
important. The fiberglass fiber diameters are approximately 10 pm, 
the G-A Kevlar 29 fiber diameters are approximately 11 pm, and the 
foam "fiber" diameters range from 31 pm through 140 pm, (75 pm was 
used in the calculations) . The calculated nonlinearity thresholds 
based on R = 1 are for the J-M 1000 fiberglass 0.11 m/s (0.1% devia- 
tion), for the G-A Kevlar 0.32 m/s (5% deviation ), and for the BAF 
0.023 m/s (Low 7% and High 12.5% deviations). 

In addition to the large "fiber” size, the low threshold of 
nonlinearity found in the BAF may also be partially explained by the 
material structure. The BAF is a partially reticulated material, 
meaning that some of the spaces between the structural members have 
thin membrane window panes. The window panes obstruct the flow and 
make the flow turbulent at a lower particle velocity than the average 
fiber size 75 pm indicates. The measured nonlinearity threshold is 
still reasonably close to the predicted value. The windows cover less 
than 25% of the foam openings. The threshold is therefore lowered only 
a relatively small amount. A material with a higher percentage of 
closed windows would be expected to have both a lower nonlinearity 
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threshold and a higher flow resistivity. This statement is supported 
by the fact that the fully reticulated Scottfelt has a high nonlinearity 
threshold and is of the same structure as the BAF but without window 
panes . 

The clc flow resistivity data presented here shows that most 
materials behave nonlinearly and illustrates the relative dc nonlinear- 
ities of the various materials. The least squares data fit is a good 
descriptor of both the linear and nonlinear regions of the flow 
resistivity. For most of the materials, Carman's analysis of the 
location of the threshold of nonlinearity is reasonable in the 
porosity region 0.8 < P< 1.0. From the previous measurement results, 
a general rule of thumb would be that Carman's modified Reynolds number 
predicts the nonlinearity threshold of fully reticulated materials to 
be in a region of between 2% and 7% of the flow resistivity deviation 
from linear. Table IV-1 lists a summary of the linear and nonlinear dc 
flow resistivity coefficients and the relative nonlinearities of the 
materials used in this study. 

2. Linear Flow Resistivity Prediction 

Hersh and Walker [ 33,34] presented two equations for use in 
predicting the linear flow resistivity of fibrous porous materials; 
see Eqs, 1-12 and 1-17, One equation can be used for calculating flow 
resistivity for flow axially along the fibers (the axial condition 
was not realized in our measurements) . Both equations can be used 
for calculating flow resistivity for flow normal to the fibers. Each 
equation has an empirical constant so that different material types 
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TABLE IV- 1 


LINEAR AND NONLINEAR DC FLOW RESISTIVITY COEFFICIENTS 
AND RELATIVE NONLINEARITIES OF BULK POROUS MATERIALS 


Material , 

Linear DC Flow 

Nonlinear DC Flow 

Relative 

Porosity P 

Resistivity 
Coefficient 
a, MKS Rayls/m 

Resistivity 
Coefficient 
n, MKS Rayls/sec 

Nonlinearii 
r i/a 

Batted 

Kevlar 

0.980 

10800. 

2370. 

0.219 

0.971 

19100. 

3860. 

0.202 

0.956 

33700. 

7930. 

0.235 

0.940 

51200. 

11400. 

0.223 

0.918 

83G00 . 

14300. 

0.172 

Scottfelt 

0.942 

13400. 

5350. 

0.399 

0.884 

49200. 

16500. 

0.335 

0.852 

66500. 

21500. 

0.323 

G-A Kevlar 

0.940 

83000. 

12000. 

0.145 

J-M 1000 
Fiberglass 

0.982 

25500. 

220. 

0.0086 

BAF 

0.970 

6390. 

27000. 

4.23 

0.971 

3900 . 

12000. 

3.08 

0.972 

8260. 

43600. 

5.28 
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may be dealt with. The two curves have been fit to measured data 

for batted Kevlar 29, as shown in Fig. IV-4. The solid curve is for 

Eq. 1-12 and the dashed curve for Eq . 1-17. The constants were chosen 

for a best Chi-squared fit to the five data points. Both curves fit 

the data points very well in the region of interest. For batted 

Kevlar 29 Hersh and Walker found that g= 0.059 fit their low frequency 

acoustic equation [ 3 3 3 ; here we find g = 0.061. Hersh and Walker 

found 1/ =0.44 [34]; here, l' = 0.445. In the first case the small 

difference in values arises from the fact that Hersh and Walker’s 

single data point was fit to an approximate acoustic equation and not 

to dc flow resistivity data. In the second case the agreement is 

excellent and supports Davies’ [24] analysis which was done with a 

variety of materials. In Fig. IV- 5 the measured dc linear flow 

resistivities for several samples of each of the other materials are 

plotted versus the porosity. The two curves, generated from Eq. 1-12 

(solid line) and Eq . 1-17 ( dashed line), show the best Chi-squared 

fits to the Scott felt data; the fits yield g= 0.0166 and ^ =1.39. 

The single data points for the other materials yield g= 0.029, and 

l / = 0.923 for the J-M fiberglass and g= 0.043, and V = 0.605 for the 
n n 

G-A Kevlar. 

The values of g and V are not the same for all materials 

n 

because of the vast differences in the material structures. The foams 
are, by no means, fibrous, the fiberglass fibers tend to bunch together 
(see Fig. III-l) and the G-A Kevlar does not have long, parallel 
strands normal to the flow. These structural effects all influence 
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FIGURE IV-4 

MEASURED DC FLOW RESISTIVITY VERSUS POROSITY FOR 
FIVE POROSITIES OF BATTED KEVLAR 29. 
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FIGURE IV-5 

MEASURED DC FLOW RESISTIVITY VERSUS POROSITY FOR GLOBE-ALBANY 
KEVLAR 29, JOHNS-MANVILLE 1000 FIBERGLASS, BLACHFORD 
ACOUSTICAL FOAM, AND THREE POROSITIES OF SCOTTFELT. 
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the way the flow is impeded by the material. Each material type must 
therefore be measured separately to determine both g and 1/ . 

In this section it has been shown that, in the region 
0.8<P< 1.0, once a dc flow resistivity, fiber diameter, and porosity 
have been determined, the change of the dc flow resistivity with 
porosity can be correlated to either Eq. 1-12 or 1-17. In the next 
section both measured and predicted dc flow resistivity values are 
used to calculate the small signal propagation of sound through the 
porous materials. 
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B. Propagation 

The results of three sets of measurements, and comparison with 
theoretical predictions, are presented in this section. Subsection 1 
is about low-intensity propagation, Subsection 2 excess attenuation and 
the approach to saturation, and Subsection 3 high-intensity propagation. 

Measurements were made on two materials, batted Kevlar 29 and 
BAF. The Kevlar is of greatest interest because of its high non- 
linearity threshold and the large range of porosities to which it can 
be compressed. The BAF is of interest because of its different struc- 
ture and low nonlinearity threshold. 

1. Propagation at Low Intensities 

Measurements of low- intensity propagation were made in a ma- 
terial filled plane wave guide, as described in Section III-D. In 
general, broadband noise in the range 0.1-10 kHz was used to determine 
the propagation parameters a and 3. The phase 3 was then used to cal- 
culate the phase speed c . Measurements made with 1, 2, and 3 kHz tone 

PH 

bursts corroberated the noise measurements. 

A highly magnified picture of batted Kevlar 29 is shown in 
Fig. III-2. The fibers are very small (-12 ym) and the material is very 
porous. The porosity range used in the propagation tests, 0. 985^_P^_0. 809 , 
was large enough to enable us to test the theory over a wide range. 

In Fig. IV-6 the predicted and measured attenuation of sound in 
batted Kevlar 29 (P = 0.980, 0.956, 0.916, 0.875, and 0.809) is plotted 
versus frequency. The horizontal bars show the range of the measured 
data points. Where no horizontal bars are drawn, either the data points 
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FIGURE IV-6 

MEASURED AND PREDICTED ATTENUATION VERSUS FREQUENCY FOR 
SOUND PROPAGATING IN FIVE POROSITIES OF BATTED KEVLAR 29. 
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coincide, or only one data point was measured. The symbols indicate the 
average value of up to seven measurements at various microphone separa- 
tions. The scatter in the measured values is due to inhomogeneity of the 
material. Since several microphone separation distances were used, the 
results therefore depended on the material homogeneity between the two 
microphones. The results of material inhomogeneity are seen later in 
the amplitude versus distance plots. 

In the theory we assumed the isothermal sound speed. In 
Appendix A this assumption is shown to be inexact. In Fig, IV-6 the 
solid line indicates the theoretical attenuation when the adiabatic 
sound speed is used in Eq. 1-4. The dashed line indicates the iso- 
thermal sound speed case. Since dc flow resistivity measurements were 
made for only five porosities (P>0.914), Eq. 1-12 (g=0.061) was used to 
calculate the unmeasured resistivities (P= 0.916, 0.875, and 0,809). 

At each porosity there are frequency ranges where the data 
and predictions agree and other ranges where they diverge. The di- 
vergence of the data above 2 kHz at porosities 0.980 and 0,916 is of 
unknown origin but was confirmed by measurements at several microphone 
positions. At the low porosities, P<0.916, there is a leveling off of 
the attenuation at low frequencies. We cannot explain this low porosity, 
low frequency attenuation. Other researchers [33,34,491 limited their 
measurements of batted Kevlar 29 to P>0.9 and did not observe this 
phenomena. 

From the arguments in Appendix A, the attenuation at the 
lowest porosity, P= 0.809 is expected to be very close to the isothermal 
(dashed) line. The data in Fig. IV-6.A show otherwise. One possible 
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explanation is that the flow resistivity o is overestimated by Eq . 1-12 
in this porosity range. We shall return to prediction of the effective 
acoustic flow resistivity later in this section. 

In Fig. IV-7 the predicted and measured phase speed of sound 
in batted Kevlar 29 (P= 0.980, 0.956, 0.916, 0.875, and 0.809) is 
plotted versus frequency. The dashed lines were calculated by using 
Eq. IV-6, The solid line were calculated by using Eq. 1-6 with the 
replacement of b Q (isothermal sound speed) with c^ (adiabatic sound 
speed). The separation between the lines indicates the limitations of 
heat transfer effects on the phase speed. Although the range bars 
appear to indicate greater error in Fig. IV-6 than in Fig. IV-7, the 
largest error, in each case, is approximately ±20%. 

We now discuss some of the limitations of these propagation 
measurement results. As the porosity is reduced the high frequency 
attenuation increases drastically, the coherence between the source and 
remote microphone signals is reduced, and the veracity of the data 
becomes questionable. The error starts in both the phase speed and 
attenuation data for p= 0.875 and f> 3 kHz. For p= 0.809 we have not 
included the f > A kHz data, because the low coherence and the large 
oscillations in the data lead to meaningless results. 

As in the low frequency attenuation results, the low frequency 
phase speed results deviate from predictions at low porosities. In the 
low frequency case, the coherence between the signals is not reduced as 
much as at high frequencies. The data deviations in Figs. IV-6 and 
IV-7 appear to be contradictory. The attenuation increase indicates 
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B. P = 0.956 AND 0.875 


FIGURE IV-7 

MEASURED AND PREDICTED PHASE SPEED VERSUS FREQUENCY 
FOR SOUND PROPAGATING IN FIVE POROSITIES OF BATTED KEVLAR 29. 
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that the phase speed should be slower, not faster than the theory pre- 
dicts. The phase speed data cannot be explained by a flanking path, 
because then the attenuation would be less than theory predicts. The 
low frequency, low porosity deviations from theory may be caused by a 
failure in the theory or in the measurement procedure. The deviation in 
question has not been resolved here. The propagation test results pre- 
sented in Subsections 2 and 3 are at 1 kHz or above; the low frequency 
deviations do not influence other results in this study. 

In Fig. IV-8 the relative sound pressure levels of three low 
Intensity tone bursts of 1, 2, and 3 kHz are plotted versus microphone 
separation distance. The material is batted Kevlar 29, of porosities 
0.980 and 0.809. The straight lines are computed by the average 
measured attenuation values (indicated by the symbols in Fig. IV-6) . Up 
to a distance of 10 cm the Fig. IV-8 data confirm the data in Fig. IV-6 
and show the materials to be homogeneous. At 20 cm the attenuation 
appears to be less than implied by Fig. IV-6. This reduced attenuation 
could be caused by one or more mechanisms, material inhomogeneity, tube 
wall vibrations, or a flanking path between the material and the tube 
wall. Although several attempts (microphone vibration Isolation and 
material repacking) were made to solve this problem, the solutions were 
unsuccessful. Nonuniform material packing is the most likely cause to 
the deviation. 

The other material tested was the Blachford Acoustical Foam 
(BAF) . The foam was cut into a long strip of square cross section and 
pulled into the tube. The predicted and measured attenuations plotted 
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FIGURE IV-8 

SMALL SIGNAL SOUND LEVELS VERSUS DISTANCE FOR THREE SINUSOIDAL 
WAVES PROPAGATING IN TWO POROSITIES OF BATTED KEVLAR 29. 



in Fig. IV-9 versus frequency show good agreement up to about 4 kHz. 

The same high frequency attenuation rise was seen in Fig. IV-6 for the 
porosity 0.980. In Fig. IV-9 the data rises to a peak and the intro- 
duction of a cross-mode propagating in the material may be causing this 
extra attenuation. 

The predicted and measured phase speeds plotted in Fig. IV-10 
show little agreement. The predicted phase speed is 40-70% too high at 
frequencies above 1 kHz. This material is a partially reticulated foam 
The structure that dc flow encounters is much different than that which 
ac flow encounters, and, for partially reticulated foam materials, the 

dc flow resistance model is inadequate for calculating 3 or C_„. 

PH 

Zwikker and Kosten £88] discuss some of the effects a partially reticu- 
lated foam will have on both dc and acoustic flow. Another theory, 
such as Zarek's [82], Beranek ? s [8], or Zwikker and Kosten T s [ 88], 
might predict the propagation properties of a partially reticulated 
foam with more accuracy than the present model. Unfortunately, in 
each of the cited theories, both the flexibility and the structure 
factor of the material must be determined, and these two factors may 
have different meaning at high-intensities. Taking these variables 
into account would also make the perturbation solution of Chapter II 
even more complicated. Thus, a partially reticulated foam is found to 
be an inapplicable case for use of the dc flow resistance theory. 
Because the relative nonlinearity of this material is so large, the 
material is used throughout the rest of this chapter as a comparator 
to the relatively linear materials. 
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FIGURE IV-9 

MEASURED AND PREDICTED ATTENUATION VERSUS FREQUENCY 
FOR SOUND PROPAGATING IN BLACHFORD ACOUSTICAL FOAM. 
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FIGURE IV-10 

MEASURED AND PREDICTED PHASE SPEED VERSUS FREQUENCY FOR 
SOUND PROPAGATING IN BLACHFORD ACOUSTICAL FOAM. 





As with the batted Kevlar 29, we want to determine how 
homogeneous the BAF is. In Fig. IV-11 the relative sound pressure 
levels of 1, 2, and 3 kHz low-intensity tone bursts in BAF are plotted 
versus microphone separation distance. In this material there is some 
deviation from the average attenuation, but it is not as large as 
found in the Kevlar material. 

In some frequency regions the measured attenuation and phase 
speed values diverge from the predicted values. If the basic model is 
valid, that is, if the attenuation is due to the material flow resisti- 
vity, agreement between theory and experiment may be obtained by using 
an acoustic flow resistivity. We measure the acoustic flow resistivity 
by using the attenuation data and Eq. 1-4 to compute a. The values of 
o are calculated by using the average of the calculations from the 
Fig. IV-6 attenuation results at 1, 2, and 3 kHz. The average values 
(#) and standard deviations are plotted versus porosity in Fig. IV-12. 
Equation 1-12 was used to generate the curve in the figure. The con- 
stant, g= 0.079, was chosen for a best Chi-squared fit to the data. For 
Kevlar flow resistivity values not measured, Eq. 1—12 was used to cal- 
culate o in the remainder of this study. 

We made the above calculations for the BAF acoustic flow 
resistivity at 1, 2, and 3 kHz and found no change between the average 
dc flow resistivity value and the specific BAF sample used in the 
propagation tests (see Fig. IV-9) . 
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FIGURE IV-11 

SMALL SIGNAL SOUND LEVELS VERSUS DISTANCE FOR THREE SINUSOIDAL 
WAVES PROPAGATING IN BLACHFORD ACOUSTICAL FOAM. 
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FIGURE IV-12 

ACOUSTIC FLOW RESISTIVITY VERSUS POROSITY FOR VALUES OF 
FLOW RESISTIVITY USED IN THE FOLLOWING PREDICTIONS OF THE 
ACOUSTICAL PROPERTIES OF BATTED KEVLAR 29. 
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2. Excess Attenuation and Saturation 


In Section I-C-l both excess attenuation and saturation 
have been shown to occur in both water [76] and air [80, 8l] for 
initially sinusoidal waves. In Section II-F-2 a model was developed 
for use in describing excess attenuation and saturation. In this sec- 
tion we present data on saturation and determine values for the nonlin- 
earity parameter T, The test method was described in Section III-D-2. 

In fitting Eq . 11-77 to the saturation data, we set the value 
of a so that Eq. 11-77 predicts a data point in the linear region and 
then set T so that Eq . 11-77 predicts the highest level data point. 
Saturation data was taken at a variety of porosities (0.980, 0.956, 
0.916, 0.895, 0.875, and 0.809), frequencies (1, 2, and 3 kHz) and 
microphone separation distances (1-20 cm). For P< 0.916 the results 
were not useable because of noise levels at the remote microphone. For 
P>_ 0.916 and several microphone separation distances the value of T 
was found to be constant with porosity changes. Approximate values of 
T at frequencies of 1, 2, and 3 kHz are listed for two materials in 
Table IV-2. A 1 kHz value of T for the BAF could not be determined from 
the experimental data because, for a 1 kHz wave in BAF, T appears to be 
dependent on sound level. This dependence is discussed more in Subsec- 
tion 3. Given these values of T, Eq . 11-77 fits almost all the data to 
within ±0.5 dB. The excellent fit is illustrated in Fig. IV-13, which 
shows the approach to saturation of three separate fundamentals for 
batted Kevlar 29 (P = 0,980). In each case, deviation from linear 
behavior begins to occur in the range of 140-150 dB source level. 
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TABLE IV-2 


EXPERIMENTALLY DETERMINED VALUES OF THE NONLINEARITY 
PARAMETER T FOR BATTED KEVLAR 29 AND BLACHFORD ACOUSTICAL FOAM 


MATERIAL 

T, 1 kHz 

T, 2 kHz 

T, 3 kHz 

Batted Kevlar 29 

388 

248 

228 

Blachf ord Acoustical 

X X 

83.6 

73.2 

Foam (BAF) 
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FIGURE IV-13 

AMPLITUDE RESPONSE CURVES FOR THREE INDIVIDUALLY 
PROPAGATING SOUNDS (1, 2, and 3 kHz FUNDAMENTAL 
COMPONENTS) IN BATTED KEVLAR (P = 0.980). 
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The corresponding particle velocity amplitude range 0. 7-2,0 m/s is 
the same range at which the dc flow resistivity begins to deviate from 
linear behavior. 

In the dc flow resistivity tests the BAF had the lowest 
nonlinearity threshold. On the basis of the results for batted Kevlar, 
one would expect BAF to have a correspondingly low threshold for extra 
attenuation. The expectation is fulfilled as Fig. IV-14, which gives 
data for two fundamental frequency components, shows. The deviation 
from linear behavior is much larger than for the Kevlar, The deviation 
starts at about 120 dB, or about 0.07 m/s in particle velocity ampli- 
tude, This amplitude also marks the onset of nonlinearity in the dc 
flow resistivity tests. 

These experiments demonstrate that the value of T depends on 
two factors, the dc nonlinearity of the material and the frequency of 
the propagating sound. In Subsection 3 we will find that, up to 172 dB, 
shocks do not form in initially sinusoidal waves propagating in bulk 
porous materials. Although shocks do not form, the equation governing 
the fundamental component amplitude attenuation is of the same form as 
for when shocks form in air. In both cases the attenuation is governed 
by an equation of the form of Eq. 11-78. 

The plot in Fig. IV-15 indicates that the apparent functional 
relation for T is inversely proportional to and (n/o)^’^. We 

find that an approximate relation is 


T = 6400(f) °- 5 (®) 0 - 4 

n 


(iv-i) 
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FIGURE IV-14 

AMPLITUDE RESPONSE CURVES FOR TWO INDIVIDUALLY 
PROPAGATING (SOUNDS 2 AND 3 kHz FUNDAMENTAL 
COMPONENTS) IN BLACHFORD ACOUSTICAL FOAM. 


138 


T 







An exact relation would be ill defined by these few tests and more 
extensive tests are called for before a good empirical relationship 
can be defined. 

Equation 11-77 and the values of T are used in the next part 
of this section to predict propagation of the fundamental component at 
high intensities. 

3. Propagation at High Intensities 

In this section we present comparisons of measured and 
predicted propagation for intense sound in porous materials. The 
propagation explored in this section is in both BAF and batted Kevlar 
29. Tone bursts of 1, 2, and 3 kHz were used to determine the material 
properties. The test method was described in Section III-D-2. 

Although propagation measurements were made on many porosities of batted 
Kevlar, only two representative data sets are presented here and com- 
pared to the theories. We find that when the perturbation solution is 
compared to the data, the solution does not adequately describe the 
propagation of intense waves In a porous material and reasons for this 
discrepancy are discussed. In addition, the impedance and saturation 
models are compared to the data. These models show the excess attenua- 
tion of the fundamental, but do not address the problem of harmonic 
component generation. 

Propagation measurements were made with batted Kevlar 29 at 
a variety of porosities (0.980, 0.936, 0.916, 0.895, 0.875, and 0.809). 
The results show general trends, which we illustrate here by presenting 
1 kHz results at only two porosities (0.980 and 0.809). The extra 
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attenuation showed that the Kevlar 29 behaves linearly for sounds 
below 140 dB. The predicted and measured sound levels of an intense 
(162.4 dB)l kHz tone burst propagating through Kevlar 29 (P = 0.980) are 
shown in Fig. IV-16. The predictions are made from a combination of 
Eqs . 11-52 and 11-61 in Eq. II-35a (Program PERT4PD) . Although the 
first four harmonics are accounted for in this solution, only three 
appear in the figure. The 4 kHz component level is less than 112 dB. 
The prediction does not show the excess attenuation that is found in 
the measured data; the computed fundamental component behaves as a 
small signal. The first two harmonic components are accounted for 
at the source microphone (0 cm) . The presence of the third harmonic 
component is not accounted for at the source microphone. The data 
points have been corrected for inhomogeneous material effects. The 
corrections are given by the differences found between data points and 
theory in Fig. TV-7. 

The data for the measured harmonic components do not follow 
the theoretical curves. The 1 kHz component attenuates in a nonlinear 
fashion. The initial attenuation is larger than linear. The attenua- 
tion approaches the linear value as the sound level decreases. The 
2 and 3 kHz components are also attenuated in a nonlinear fashion. 

The 2 kHz component has a rapid initial attenuation which is reduced 
with distance. The 3 kHz component increases in level at the 2 and 
5 cm positions, but is suddenly attenuated between 5 and 11 cm. These 
results indicate interaction between the components which is unex- 
plained by the theory. The discrepancies may be caused by one or more 
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FIGURE IV-16 

MEASURED AND PREDICTED (PERT4PD) HARMONIC COMPONENTS SOUND 
LEVELS VERSUS DISTANCE FOR AN INTENSE 1 kHz SOUND (162.4 dB) 
PROPAGATING IN BATTED KEVLAR 29 (P = 0.980). 



factors. First, the presence and relative phases of the higher fre- 
quency harmonic components may influence the lower frequency components 
more than expected. Second, the perturbation may have to be performed to 
higher orders to show enough excess attenuation. Third, the definition 
of nonlinear behavior may be ill defined by simply using the dc flow 
resistivity test results; a better definition may be needed. Finally, 
perturbation may be inappropriate for determining the propagation of 
very intense sound in a porous material. We found in Section IV-B-1 
that for small signals the perturbation theory and measurement results 
were in agreement. Until harmonic generation by and excess attenuation 
of the fundamental is encountered, the perturbation theory works well in 
predicting the attenuation (below about 140 dB) . 

We now consider how well the perturbation theory predicts the 
attenuation of sound for materials of low porosities. In Table IV -1 we 
found that the relative nonlinearity of the material n / cr is independent 
of the porosity. The resistivity and attenuation are dependent on the 
porosity. At low porosities G-*0 (Eq . 11-81) and the nonlinearity is 
overshadowed by the large attenuation. In Fig. IV- 17 the perturbation 
solution and data are presented for batted Kevlar 29 P= 0.809. Here 
there is very little excess attenuation; theory and experiment agree. 

The 2 kHz component shows linear attenuation between the two data 
points. A small rise is found in the predicted 2 kHz level. The dis- 
crepancy between data and theory here may be caused by lack of initial 
phase information. Changing the phase in the computer program changed 
the generation and attenuation of the 2 kHz component. The phase in 
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FIGURE IV-17 

MEASURED AND PREDICTED (PERT4PD) HARMONIC COMPONENT SOUND 
LEVELS VERSUS DISTANCE FOR INTENSE 1 kHz SOUND (161.8 dB) 
PROPAGATING IN BATTED KEVLAR 29 (P = 0.809.). 




Fig. IV-17 was set to zero. For porosities between the two examples, 
we found a transition from some excess attenuation at high porosities 
to little excess attenuation at low porosities. 

We now use test results to verify the two attenuation models 
presented in Section II-F, We use the values of the nonlinearity param- 
eter T from Section IV-B-2 and the values of the relative nonlinearity 
n/a from Section IV-A-1. The acoustically determined flow resistivity 
values are taken from Section IV-B-1. We first consider the 1 kHz 
Kevlar data and then the 1 kHz BAF data. We find that the measure- 
ments and predictions agree for the Kevlar, but there are some discre- 
pancies for the BAF. 

In Fig. IV-18 the measured and predicted attenuations of a 
1 kHz, 162.4 dB, wave are plotted versus distance. This is the 
same data as used in Fig. IV-16. The solid line indicates the small- 
signal attenuation. The dashed line is a plot of Eq . 11-77, the 
excess attenuation model and the dotted line is a plot of the result 
of using Eq. 11-70, the impedance model. The excess attenuation model 
fits the data the best. The impedance model predicts too much excess 
attenuation, a fact that indicates that, for the impedance approxima- 
tion the values of n/a derived in Section IV-A-1 are too large for the 
approximation. The need to reduce the values of n/a is confirmed by the 
other comparisons of data with theory, including those of impedance. 
Since the predictions based on Eq. 11-79 show good agreement with the 
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FIGURE IV-18 

MEASURED AND PREDICTED FUNDAMENTAL COMPONENT SOUND LEVELS 
VERSUS DISTANCE FOR AN INTENSE 1 kHz SOUND (162.4 dB) PROPAGATING 

BATTED KEVLAR 29 (P = 0.980). 




data, Eq . 11-79 will be used in the remainder of this section. 
Estimates for the reduction of r\/o are made in Section IV-C. 

In Fig. IV-19 we continue with the very porous batted Kevlar 
29 (P= 0.980) and increase the sound level at the source microphone 
to 172.3 dB, again for a fundamental freuqency of 1 kHz. The conclu- 
sions are the same: predictions based on the use of Eq . 11-79 are 

in excellent agreement with the data over the region measured. 

A similar measurement in low porosity (P = 0.809) batted 
Kevlar 29 is shown in Fig, IV-20. In this case the material nonlin- 
earity is not as important as the linear attenuation. Even so, the 
agreement between data and predictions based on Eq . 11-79 continues 
to be excellent. Moreover, the earlier conclusion that T for a given 
material is independent of porosity is now seen to hold for porosities 
as low as P = 0 . 809 . 

The same sort of tests were done with BAF. It will be 
recalled that a 1 kHz value for T could not be determined from the 
saturation measurements, even though the 2 and 3 kHz values could be 
determined. From Eq . 11-79 we estimate 7= 115. In the BAF tests 

for 162.1 dB presented in Fig, IV-21, we find T 1 =115 to be a good 
estimate. In Fig. IV-22 the initial sound level is 173.1 dB. We 
find that the model and data do not agree as well. The estimate for 
the BAF value T=1I5 appears to be low. If T is increased to fit the 
two cases shown here, the results are T=143 for Fig. IV-21 and T = 200 
for Fig. IV-22. In other words, T seems to have an amplitude dependence 
that was not seen for Kevlar 29. The dependence on amplitude may be 
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FIGURE IV-19 

MEASURED AND PREDICTED FUNDAMENTAL COMPONENT SOUND LEVELS 
VERSUS DISTANCE FOR AN INTENSE 1 kHz SOUND (172.3 dB) PROPAGATING 

IN BATTED KEVLAR 29 (P = 0.980). 
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caused by the complicated material structure which influences the 
acoustic flow in a different manner as the intensity increases. For 
the 2 and 3 kHz tests amplitude dependence of T was not encountered. 

We have seen here that, although the second-order perturba- 
tion solution does not work for intense sound propagation in a porous 
material, we have two empirical models that do work. Future work in 
theoretically defining the parameters T and n/a would be an important 
step in being able to predict the propagation of intense sound in 
porous materials. 

In Section IV-C the acoustic impedance theory and data are 
compared. Acoustic values for n/a are found for use in the impedance 
relations . 
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C. Specific Normal Acoustic Impedance 

The acoustic impedance measurement systems used in this 
section are described in Section III-E. The measurement systems con- 
sist of two low-intensity sound standing wave tubes and one high- 
intensity sound traveling wave tube. Measurements were done in these 
tubes for a wide range of frequencies, materials, and sound levels. 

All the results cannot reasonably be reported here. Tabel IV-3 lists 
the range of measurements made on the various materials. Low-frequency 
standing wave tube measurements were not made on the batted Kevlar 29 
because of our inability to hold the material in place at various 
porosities and keep the material flush with the tube walls. The 
low-intensity measurements are presented and discussed first, then the 
high-intensity measurements, 

1. Impedance at Low Intensities 

Most porous materials are of finite length and the measured 
impedance is affected by the material impedance and the sample length. 
The impedance data presented here were obtained from both the standing 
and traveling wave tubes. In the standing wave tube the sample lengths 
were of the order of 7.6 cm. In the traveling wave tube the sample 
was long enough to qualify as semi-infinite (aL>>l). We use the 
traveling wave tube to measure the characteristic impedance of the 
material. The test results from the traveling wave tube are presented 
first. The flow resistivity values used in the calculations are taken 
from both the dc flow resistivities, described in Section A of this 
chapter, and the acoustically determined flow resistivities, described 
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TABLE IV-3 

TABULATION OF MEASUREMENTS MADE IN THE THREE 
IMPEDANCE TUBES ON THE VARIOUS MATERIALS 


MATERIAL 


BATTED KEVLAR 29, 


POROSITY 

TRAVELING 
WAVE IMPEDANCE 
TUBE 

-(FREQUENCY, kHz) 
-(SOUND LEVEL, dB) 


0.985 

1-4 

119-164 


0.971 

1-4 

121-167 


0.958 

1-4 

120-167 


0.936 

1-4 

121-167 


0.914 0.895 0, 

1-4 1-4 ! 

121-166 121-164 121 


LOW FREQUENCY 
STANDING WAVE 
IMPEDANCE 
TUBE 

-(FREQUENCY, kHz) 


HIGH FREQUENCY 
STANDING WAVE 

IMPEDANCE TUBE " 0.8-3. 5 - 0.8-3. 5 - 0.8-3. 5 0. 8-3.5 

-(FREQUENCY, kHz) 


MATERIAL 

JOHNS-MANVILLE 
1000 FIBERGLASS 

GLOBE -ALBANY 
KEVLAR 

900-Z-2 

SCOTTFELT 

900-Z-4 

900-Z-6 

BLACHFORD 

ACOUSTICAL 

FOAM 

POROSITY 

0.982 

0.940 

0.942 

0.884 

0.852 

0.971 

TRAVELING 
WAVE IMPEDANCE 
TUBE 

-(FREQUENCY, kHz) 
-(SOUND LEVEL, dB) 

- 

1-4 

120-165 

1-4 

120-165 

1-4 

120-165 

- 

1-4 

120-165 

LOW-FREQUENCY 
STANDING WAVE 
IMPEDANCE TUBE 
-(FREQUENCY, kHz) 

0.1-1. 3 

- 

- 

- 

- 

0. 1-1.3 

HIGH-FREQUENCY 
STANDING WAVE 
IMPEDANCE TUBE 

0.8-3. 5 

0.8-3. 5 

0.8-3. 5 

0.8-3. 5 

0.B-3.5 

0.8-3. 5 


-(FREQUENCY, kHz) 


875 

-4 

-165 
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in Section B of this chapter. The characteristic impedance data il- 
lustrate how well the low-intensity sound model predicts the impedance 
for a semi-infinite material. The measurements are combined with the 
theory of Chapters I and II to yield a prediction for the impedance of 
the finite materials. 

For a finite length porous material backed by an infinite im- 
pedance a simple relation [88] was found. This relation, given by 
Eq. 1-10, links the propagation parameters a and 8, the semi-infinite 
material impedance, and the finite length material impedance W(L) 
together. Equation 1-10 can be rewritten in terms of the dimensionless 
variable F ^ as 

/ jcur i L \ 

Z 1 (L) = T 1 coth (— ^-1 (IV-2) 

We note that, when aL>> 1, - F^. As noted in Chapter I, the measure- 

i 

ment of the small signal impedance of a semi-infinite material is 
equivalent to measuring the propagation parameters. 

In Fig. IV-23 the predicted (Eq. 11-62) and measured (travel- 
ing wave impedance tube) impedance of batted Kevlar 29 (P - 0.985, 0.956, 

0.914, and 0.875) are plotted versus frequency. For the higher poros- I 

ities the data and theory are in good agreement. At P= 0.875, however, 

as in the attenuation and phase speed tests, the predicted values are 

larger than the measured values. The acoustically measured resistivity 

values are used in the calculations. If the dc flow resistivity values ^ 

are used, the predicted curves for the impedance are higher than the curves 

in Fig. IV-23. At the higher porosities (P > 0.875) the small discrepancies 
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FIGURE IV-23 

MEASURED AND PREDICTED IMPEDANCE VERSUS FREQUENCY FOR 
FOUR POROSITIES OF SEMI-INFINITE BATTED KEVLAR 29. 



between the measured and predicted values can be accounted for by the small 
inaccuracy in measuring the phase (time delay) in the tests. 

Samples of 7.7 cm thick batted Kevlar 29 were measured in 
the standing wave tube. The data are compared with predictions 
based on Eq. IV-2 (acoustic resistivity values) in Fig. IV-24 for 
porosities 0.980 and 0.918. Except for the P = 0.918 resistance data, 
the agreement is seen to be good. When the dc flow resistivity values 
are used in the predictions, the differences are barely discernable 
from the curves in Fig. IV-24. It would appear that for aL < 1, that 
the differences between how the flow resistivites are determined 
becomes unimportant and either flow resistivity may be used in the 
calculations . 

We stated that when aL > 1, the standing and traveling wave 
impedance tube measurements should yield similar results, A comparison 
of Figs. IV-23 and IV-24 indicates that the test results should be 
similar for P =0.918, 0.914, where aL > 1 and dissimilar for P= 0.980, 

0.985, where aL < 1. The theoretical curves indicate the above state- 
ment to be fairly reliable. The data for P= 0.985, 0.980 disagree in 
the manner indicated by the theory. The reactance data for P= 0.918, 

0.914 agree, but the resistance data disagree. The resistances 

(P = 0.918, 0.914) measured in the standing wave tube are between 15 and 35% 

higher than both the theory and the resistances measured in the traveling 

4 ! ‘ 

wave tube. As we shall see, the disagreement in the measured resistances 
for P= 0.918 and 0.914 appears to be more in this specific data than in 
the procedure, i.e., we find more consistant agreement in other results. 
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FIGURE IV-24 

MEASURED AND PREDICTED IMPEDANCES VERSUS FREQUENCY FOR 
TWO POROSITIES OF 7.7 cm THICK BATTED KEVLAR 29. 



In Fig. IV-25 the predicted and measured impedances of 


needled and felted Kevlar 29 (P= 0.940) for semi-infinite and 6.2 cm 
samples are plotted versus frequency. No propagation tests were made 
on this material; the dc flow resistivity was used in the predictions. 
The two sets of high-frequency data are in good agreement. The agree- 
ment between theory and data is good at high frequencies but not at low 
frequencies. The discrepancy at low frequency is attributed to un- 
certainty in locating the impedance tube microphone, a problem discussed 
in Chapter III. 

In Fig. IV-26 the impedance of a 6.2 cm thick sample of 
Johns-Manville 1000 fiberglass (P = 0.982) is plotted versus frequency. 
The agreement between theory (dc flow resistivity) and data is excel- 
lent, significant deviation occurs at only the lowest frequencies. 

The next material is Scottfelt 900-Z-2 (P = 0.942). The 
predicted (dc flow resistivity) and measured impedances of the Scott- 
felt for semi-infinite and 7.5 cm samples are plotted versus frequency 
in Fig. IV-27. The measured resistances are between 0 and 50% higher 
than predicted. The measured reactances are in good agreement with 
the predictions. In both cases the measured results are in excellent 
agreement. The agreement between theory and experiment for other 
Scottfelt samples is similar. The dc flow resistivity was used in 
the theoretical predictions. Although agreement between theory and 
data is not quite as good as for the fibrous materials, it is still 
fairly good. We therefore conclude that the behavior of the fully 
reticulated foam is similar to that of a fibrous material. One reason 
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FIGURE IV-26 

MEASURED AND PREDICTED IMPEDANCES VERSUS FREQUENCY FOR 
6.2 cm THICK JOHNS-MANVILLE 1000 FIBERGLASS. 
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this material was not used in the propagation parameter tests is that 
many pieces were not uniform in porosity throughout their thickness. 

The pieces were, in many cases, less porous at the surface than in the 
center. Although the samples chosen for this test were apparently 
uniform, invisible nonuniformity may have been one of the causes of 
the higher measured than predicted acoustic resistance. 

In Fig. IV-28 the predicted and measured impedances of BAF 
(P= 0.971) for semi-infinite and 7,6 cm samples are plotted versus fre- 
quency. As expected from the propagation test results, the agreement 
between the theory and data is not good. The disagreement is much 
greater than for Scottfelt, but the trends are similar: The resistive 

part is much higher than predicted, whereas the reactive part is close 
to the predicted curve. 

By comparing the data and theoretical results in Figs. IV-23 
through IV-28, we find that fibrous and fully reticulated foam ma- 
terials have similar acoustical properties and may be treated as rigid 
materials in the small signal case. On the other hand, the partially 
reticulated foam (BAF) has a measured acoustic resistance that is in- 
explicably high. As noted in Section IV-B above, one of the more 
general theories (see Chapter I) may provide a better explanation of 
the behavior of a partially reticulated foam. 

2. Impedance at High Intensities 

Proceeding with the ongoing quest for the definition of a 
useful nonlinearity parameter for bulk porous materials leads to the 
ad hoc nonlinear impedance model presented in Section II-E. The 
model consists of substituting the nonlinear dc flow resistivity 
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FIGURE IV-28 

MEASURED AND PREDICTED IMPEDANCES VERSUS FREQUENCY FOR 
7.6 cm THICK AND SEMI-INFINITE BLACHFORD ACOUSTICAL FOAM. 




a^ + £ri^/Pv^ for the linear flow resistivity in the definition of 
High-intensity measurements were not made in the standing wave 
tubes and only the test results for semi-infinite materials are 
presented. The materials are the same as those reported in part 1. 

In making many comparisons between the measured and pre- 
dicted nonlinear impedance behavior, we find that use of the dc flow 
resistivity in the theoretical calculations leads to an incorrect 
prediction of the material nonlinearity. (In the cases where acoustic 
resistivity values were used in the predictions, the same ratio of 
rj/c as found in Table IV-1 was employed to calculate n.) For the 
following comparisons the values of n are fit to the data. 

In Fig. IV-29 the impedance at 1 kHz of several porosities 
(P= 0.985, 0,958, and 0.895) is plotted versus sound pressure level. 
The relative nonlinearities used in this figure are n/a=0.24, 0.12, 
and 0.06, whereas in Table IV-1 we found n/a-0.2. As with the ex- 
cess attenuation parameter T and its relationship to the attenuation 
a, we find that the importance of the nonlinearity apparently de- 
creases with increased attenuation. 

Neither the resistive nor reactive parts of the impedance 
vary in magnitude to any great extent. In each of the three cases the 
measured reactive part appears to be slightly more linear in behavior 
than the resistive part. This behavior is found to exist at all 
porosities and frequencies in this material. The asymmetric behavior 
is illustrated in Fig. IV-30 where the frequency dependence of the 
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FIGURE IV-29 

MEASURED AND PREDICTED IMPEDANCES AT 1 kHz VERSUS SOUND PRESSURE 
LEVEL FOR THREE POROSITIES OF SEMI-INFINITE BATTED KEVLAR 29. 



FIGURE IV-30 

MEASURED AND PREDICTED IMPEDANCES AT TWO SOUND LEVELS VERSUS 
FREQUENCY FOR SEMI-INFINITE BATTED KEVLAR 29 (P - 0.936). 






batted Kevlar (P = 0.936, ri/a = 0.12) impedance is plotted at two sound 
levels. The resistive part of the impedance is accurately predicted 
for all frequencies at both low and high intensities. On the other 
hand, the reactive part is well predicted at low intensity, but not 
at high intensity. Constancy in the reactive part of the impedance 
might be caused by fiber motion or by near independence of the 
reactance on particle velocity. 

In Fig. IV-31 the predicted and measured 1 kHz impedances 
of the needled and felted (G-A) Kevlar 29 (P= 0.940) and the Scottfelt 
(P= 0,942) are plotted versus sound pressure level. The porosities of 
the materials in this figure are similar to the material in Fig. IV-30, 
but the material structures are quite different. The relative non- 
linearities fit to the data are ri/a = 0.08 for the G-A Kevlar and 
T)/o = 0.44 for the Scottfelt. 

The last material considered is the BAF , The dc flow 
resistivity tests indicated that this material has very nonlinear 
behavior. This fact is supported in the plot of the 1 kHz impedance 
versus sound level in Fig. IV-32. The relative nonlinearity r \ /a = 1.2 
has been fit to the data. The relative nonlinearity is much less than 
that listed in Table IV-1 (n/o~4). As opposed to the predictions for 
the other materials, the trends of both the resistive and reactive 
parts of the impedance are predicted by the nonlinear impedance 
equations. Except for the last data point, the resistive part is 
slightly more nonlinear than the reactive. A conclusion we draw from 
this figure is that the impedance model predicts changing reactance 
for all materials. In the cases of fibrous or fully reticulated 
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FIGURE IV-31 

MEASURED AND PREDICTED IMPEDANCES AT 1 kHz VERUS SOUND PRESSURE 
LEVEL FOR GLOBE-ALBANY KEVLAR 29 AND SCOTTFELT 900 Z-2. 
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FIGURE IV-32 

MEASURED AND PREDICTED IMPEDANCES AT 1 kHz VERSUS SOUND 
PRESSURE LEVEL FOR SEMI-INFINITE BLACHFORD ACOUSTICAL FOAM. 




foam materials the threshold of nonlinear reactive behavior may be 
higher than anticipated by this model. 

In Fig. IV-33 the frequency dependence of -BAF at two sound 
levels is shown. We again see that the behavior is quite different 
from the batted Kevlar (Fig. IV-30) and that while the reactive part 
is accurately predicted at both sound levels, the resistive part is 
not. 

In this ad hoc model, as opposed to the other high intensity 
sound models presented, the nonlinear effects are apparently dependent 
on the porosity (and, in turn, the attenuation). The nonlinearity de- 
creases in importance as the porosity decreases, but this has not been 
accounted for in the model. A specific relationship between porosity 
and relative nonlinearity has not been drawn because the relationship 
appears to differ with the material. 







D. Summary 

In this chapter the theoretical predictions have been com- 
pared with data from experiments for flow resistivity, propagation 
parameters, saturation effects, and specific normal acoustic impedance. 
Where possible, analysis has been presented with respect to frequency, 
sound pressure level, porosity, and material type. The agreement and 
disagreement between the data and theory have been discussed. The 
major points of the analysis are summarized in this section. 

Two conclusions from the dc flow resistivity measurements 
(Section A) are drawn. First, each material exhibits nonlinear be- 
havior at some threshold particle velocity. The onset of nonlinearity 
appears to be related to the material structure. For most fibrous 
materials and fully reticulated foams the modified Reynolds number can 
be used to estimate the nonlinearity threshold. Second, the dc linear 
flow resistivity was found to be a function of the porosity and the 
material structure. The exact relationship is unknown, but two pre- 
viously derived functions [33,34] Eqs, 1-12 and 1-17 fit the data for 
two materials very well. Although only single porosities of some 
materials were measured, the functions should be good descriptors of 
the porosity-flow resistivity relationship. 

The propagation measurements (Section B) yield important in- 
formation about the definition of acoustic material nonlinearity. 
Analysis of the small-signal measurements (Part 1) showed that use of 
the dc flow resistivity data leads to a small, systematic error in the 
prediction of the acoustic propagation parameters of the porous ma- 
terials. In the case of the batted Kevlar 29 material, the empirical 


173 



constant g in Eq, I 12 was changed to reflect a more realistic value 
for the material properties. No changes were made in the empirical 
constants of the other materials. 

In Part 2 saturation effects in bulk porous materials were 
studied. It was determined that, although shocks never form, at high 
intensities the energy loss by the fundamental and the approach to 
saturation are governed by the same mathematical relationship in 
^‘3* 11-76 as when shocks do form in waves propagating in fluids. A 
Gol dberg-type number is found in Eq . 11-79. From comparison with data 
we define values for a porous material nonlinearity parameter T which 
is part of the Gol' dberg-type number. The nonlinearity parameter was 
found to depend on frequency and relative nonlinearity (material type), 
but, at least for batted Kevlar 29, not on porosity. High-intensity 
propagation measurements were reported in Part 3 and compared with the 
results of the perturbation solution of the conservation equations. At 
low intensities the propagation of the sinusoidal waves was found to be 
affected by material inhomogeneities. At high intensities the in- 
homogeneities affected the results and could be corrected for. At the 
higher porosities of Kevlar 29 the second-order perturbation solution 
results were found to diverge from the measured data. The theory 
predicted too much harmonic generation and not enough excess attenua- 
tion of the fundamental. As the porosity was decreased the theory was 
found to agree with the data at higher sound levels. The energy lost 
by the fundamental to excess attenuation was, apparently, not all 
shifted to the higher harmonics, as it is in fluids. It appears that 
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shocks do not form in the porous material because the nonlinearity is 
a turbulence related phenomena and the transference of energy out of 
the fundamental may be done on a broadband as well as a pure tone basis. 
There is no experimental basis for this statement, but it is offered as 
a possible explanation to where the energy goes. 

Section C (impedance measurements) contains important infor- 
mation for the definition of the material impedance and support in- 
formation on the propagation parameters. In Part 1 the low intensity 
theoretical and measured impedance for semi-infinite and finite ma- 
terials were found to agree. Since the propagation parameters and the 
impedance of semi-infinite materials are directly related to each 
other, the agreement between data and theory supports the calculations 
and measurements for the propagation parameters. 

In Part 2 the ad hoc impedance model was compared to the im- 
pedances of several materials at different frequencies and sound levels. 
The definition of the material nonlinearity through the use of this 
model acts in opposition to the dc flow resistivity model as well as 
the saturation model. The nonlinearity defined by the nonlinear impe- 
dance model is independent of frequency, but depends upon the porosity 
(attenuation) of the material. This result is not surprising because 
there is no actual theoretical basis for combining second order effects 
in a first order impedance equation. As a rough approximation, the 
approach works. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

The propagation of intense sound in bulk porous materials 
has been studied both theoretically and experimentally. In the theo- 
retical analysis, rigid material structure has been assumed and the 
viscous loss has been defined by the measured dc flow resistivity. In 
the linear analysis region these analytical assumptions have been 
successfully used by others [23,32,33,34,45] to explain their data. 

We have also experimentally verified the theory here. 

A brief, small signal analysis of oscillatory heat transfer 
effects in fibrous porous materials is presented in Appendix A. Heat 
transfer effects modify the compressibility of the gas and, thus, the 
reference sound speed. The resultant reference sound speed is a 
function of fiber size, porosity, and frequency. Because the heat con- 
ductivity and heat capacity of porous materials are much larger than 
that of the air, the type of material used does not affect the analy- 
sis. The heat transfer analysis results are not used directly in the 
analysis presented in the main body of this study because the effects 
are small relative to the viscous effects. The heat transfer effect 
limits (isothermal and adiabatic sound speeds) are used to indicate 
the range of modifications to the viscous theory. 

Standard nonlinear acoustic theory has been used as a basis 
for understanding the nonlinear effects on intense sound propagating in 
porous materials. Unfortunately, nonlinear acoustic theory of fluids 
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is not directly applicable to intense sound propagation in a porous 
material because of the high attenuation and severe dispersion found 
in these materials. For sounds propagating in fluids, viscous attenua- 
tion and dispersion effects are usually ignored in the pre-shock re- 
gion. In bulk porous materials the attenuation and dispersion are so 
large that shocks never form. Despite the lack of shock formation, 
some of the effects of shock formation on the fundamental component, 
excess attenuation and saturation, are found to occur as if shocks had 
actually formed. 

The experimentally determined approach to saturation is 
useful in describing how much of the acoustic energy of intense sound 
is propagated, or not propagated, in the porous material. The modeling 
of saturation effects was based on knowledge of how the energy is lost 
from the fundamental in the sawtooth region of an intense wave in a 
fluid [76,80,81]. For the high intensity energy loss an empirical 
parameter T has been defined and found to depend upon frequency f and 
the relative nonlinearity r\/a of the material. Because of the limited 
data, the exact relationship T(f,n/c) is unknown, but we defined an 
approximation. It is possible that T is a function of porosity. The 
excess attenuation and saturation experiments might be used to define 
the acoustical nonlinearity of the material. Possibly, a relation- 
ship between the saturation defined acoustic nonlinearity and the dc 
flow resistivity nonlinearity could be determined. There is probably 
a definable function between the two nonlinearities for fibrous mate- 
rials or fully reticulated foams, but it may not be reasonable to 
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expect such a function for partially reticulated foams. In most 
cases we find that the fit of the excess attenuation model to the 
data to be excellent. 

The perturbation solution of the conservation equations was 
derived in an attempt to explain excess attenuation of the fundamental 
and harmonic component generation in the measurement results. The 
model does not work well at high Intensities, porosities, and 
nonlinearities. As the porosity is lowered the agreement between 
prediction and measurement results is improved. 

Success of the excess attenuation model and failure of the 
perturbation theory have led us to re-examine how the excess attenua- 
tion energy is lost from the fundamental component. It appears 
that the fundamental component energy is lost not only to the harmo- 
nic components, but, because the nonlinearity is a turbulence process, 
to broadband random noise. This broadband noise level can be as 
much as 60-70 dB down from the fundamental component level, thus 
making it difficult to measure and not noticeable. 

Impedance measurements and predictions for many materials 
have been presented to illustrate that the theory is applicable to 
most finite and semi-infinite materials at low intensities. These 
results are also presented In support of the propagation parameter 
measurements, since when aL>>l, + and the same results are 
obtained for each measurement. 

A new experimental technique was developed for use in measur- 


ing the Impedance of semi-infinite porous materials at high intensities. 



The technique has some advantages over techniques developed by others. 
The system is simple to make and use and measurements may be made up to 
sound levels of 167 dB. As with some other measurement techniques, 
shock formation limits the sound level to which this method is appli- 
cable. 

High-intensity sound measurements were done to evaluate how 
well an ad hoc model predicts the effect of sound level on impedance 
of a semi-infinite material. For the ad hoc model to be useful, the 
relative nonlinearity of the material has to be reduced with reduced 
porosity in order to predict the measured values. Because of the 
simplicity of the saturation and impedance models, relationships be- 
tween the three experimentally determined nonlinearities would be 
useful in acoustical analysis. The determination of these relation- 
ships appears to be a fruitful area for future research. 

The nonlinear effects described in this study are only 
applicable to materials exposed to very intense sound fields; even 
the most nonlinearly behaving materials do not exhibit nonlinear 
behavior below 120 dB. The most important area for use of these re- 
sults is to bulk material use in the sound reduction of jet engine 
noise. Materials such as Kevlar 29 have recently been used inside the 
Helmholtz resonators that act to absorb sound before it is radiated 
out the engine intake. The additional treatment is, apparently, 
successful in reducing the radiated noise [92]. There may be a com- 
bination of effects caused by the porous material that offsets the 



nonlinear behavior of the resonators and increases their effectiveness 


relative to resonators without porous materials. For example, the 
porous material impedance changes may offset the resonator impedance 
changes with sound level. Or, saturation effects may increase the 
sound absorption of a sound once it has entered the resonator, not 
letting as much sound out as without the porous material. 

This study has been successful in describing how sinusoidal 
(or almost sinusoidal) waves interact with a bulk porous material. 

The results of the study should also be useful in laying the ground- 
work for future studies in this area. Important topics for future 
research would be (1) an investigation of the relationships between the 
acoustic nonlinearities as defined in the various theories and experi- 
ments in this study and (2) an investigation of the interaction of the 
acoustic field and a resonator filled with a fibrous material. 
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APPENDIX A 


A Theoretical Description of Oscillatory Heat Transfer Effects 
on the Propagation of Sound in a Fibrous Porous Material 

In this appendix a mathematical model of heat transfer effects 
on sound propagation in air saturated fibrous material is derived. The 
resulting model is good for low- intensity sounds propagating in the air 
of the porous material. The conductive heat transfer properties of both 
the air and the fibers are accounted for. In the final analysis, the 
specific heat transfer properties of the fibers are found to have a 
minimal effect on the sound propagation and may be neglected. The heat 
transfer primarily affects the sound speed. Its influence on the atten- 
uation is small relative to that of viscosity and is ignored. 

The derivation that follows is for cylindrical fibers, and 
cylindrical coordinates are used. The heat transfer is assumed to be 
due to conduction only, i.e., convective heat transfer and acoustic 
streaming (mixing) at high intensities are not accounted for. A primary 
effect of convection and streaming would be to raise the heat transfer 
rate and thus increase the tendency of the propagation to be isothermal. 

Two steps are used in determining the oscillatory heat trans- 
fer in fibrous porous materials. First, the general problem of oscilla- 
tory heat transfer between the medium and the cylindrical fiber is 
determined. Two different boundary conditions at the fiber surface are 
considered, and the effects of each on the sound propagation are dis- 
cussed. Second, the relation describing the effect of the heat transfer 
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on the phase speed is derived. This relation is used in Eq . 1-6 for 
comparison to the phase speed of Kevlar 29 (P = 0.980). 

I. Oscillatory Heat Transfer in Cylindrical Coordinates 

In this section the oscillatory heat transfer between a single 
fiber and the infinite material in which it is immersed is considered. 
The general heat transfer equation is presented and solved for two sets 
of boundary conditions. The two solutions are compared. 

A. Equation and General Solution 

The governing equation for radial heat conduction in cylin- 
drical coordinates is [22] 



+ 1 
r 3r/ 


J_ 3T 

t* T 8t 


(A-l) 


where a^, = K/pC is the thermal diffusivity, p is the density, C is the 
specific heat, k is the thermal conductivity, r is the radial distance 
from the fiber center, T is the temperature, t is the time, and A is an 
arbitrary driving function. The driving function A may be a function of 
the coordinate system and time. The relation kA is the rate of heat 
production per unit volume [22]. 

If the driving function A is sinusoidal in time, i.e., 

A = h Q e iuit , (A-2) 

then the steady state solution of Eq. A-l is [22,52,61] 


T(r,t) 



+ DK 

o 






+ T 
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(A-3) 
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where J is the Bessel function of the first kind of order zero, K is 
o o 

the modified Bessel function of the second kind of order zero, is 
the reference (ambient) temperature, and £ = /co/a , In this presenta- 
tion we assume that A is constant in the region of interest, either 

o 

internal or external to the fiber, Bessel functions with imaginary 
arguments can be written in terms of Kelvin functions with real argu- 
ments. The relations are [22,52,6l] 

( z J 3j/2 ) = k er v z + j bei^z (A-4) 

and 

K v ^zj 1/2 j = j V ^ker^z + j kei^z) . (A-5) 

Alternately, the Kelvin functions can be written in polar form as 


(z) e 


j 0 v (z) 


yjber^z + bei^z ^cos0 y (z) + j sinS^Cz 


(A-6a) 


and 


N v (z)e 


j<t> v (z) = 


yjker^z + kei 2 z ^cos^^(z) + j sin<j> v (z)^ , (A-6b) 


where 


bei z 

0 (z) = arctan - — — (A-7a) 

v ber z 

v 


18 3 


and 


kei z 

♦v (z) = arctan • (A-7b) 

V 

The equations above are valid in most regions of the coor- 
dinate system. The K^zj term becomes infinite at z=0 and the 

/ 3/2 \ 

J^zj j term becomes infinite at z=°°. In what follows, the subscript 
I denotes the region internal to the fiber and the subscript E the 
surrounding material. The constants B and D are determined by the 
boundary conditions, which are discussed next. 

Table A-l is a listing of the thermal properties of various 
materials. Several of these values are used in following examples. 

B. Boundary Conditions 

Two sets of boundary conditions at the fiber surface are con- 
sidered here. In one case the fiber surface is assumed to be isothermal, 
in the other case non-isothermal. But first, general boundary con- 
ditions, that is, conditions at r=0 and 00 must be stated. 

1 , General 

The boundary condition at the fiber center is that the tem- 
perature is finite, or 

= 0 . (A-8) 

r-0 


3T E (r,t) 

9r 


The boundary condition at an infinite distance from the fiber is 
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TABLE A-l 


THERMAL PROPERTIES OF SEVERAL MATERIALS 


MATERIAL 

K 

THERMAL 

CONDUCTIVITY 

W/m°K 

p 

DENSITY 

kg/m^ 

SPEC?FIC 

HEAT 

J/kg°K 

THERMAL 

DIFFUSIVITY 

m^/s 

f =1000 1 
a= 6 - 10" 6 

COPPER [50] 

390 

8900 

385 

-4 

1 . 1*10 

0.045 

STEEL [50] 

45 

7700 

450 

1.3* 10 ~ 5 

0.13 

AIR (300°K) 
[28] 

0.026 

1.2 

1006 

-9 

7.5-10 

0.10 

GLASS [50] 

2.5 

2600 

775 

1 . 2 ’ 10~ 6 

0.42 

KEVLAR [50] 

0.7 

1445 

1200 

4.0-10 " 7 

0.75 

NYLON [50] 

0.2 

1100 

1700 

i. 1 ’ 10 -7 

1.5 

TEFLON [50] 

0.24 

2200 

1050 

l.o-io -7 

1.5 

POLYURETHANE 

0.3 

1150 

1000 

2 . 6 - 10" 7 

0.93 


[50] 
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(A-9) 


9T £ (r,t) 

9r 


= 0 

r =°° 


2. Isothermal Surface 

The mathematical analysis is simplest if the surface tempera- 
ture is constant, i.e., isothermal, 

T(a,t) = T , (A- 10) 

o 

where a is the fiber radius and is the reference (ambient) tempera- 
ture, In this case the external and internal temperature fields are 

i 

decoupled. As will be seen, this condition closely approximates the 
actual surface boundary condition of the fiber In air. The isothermal 
boundary condition is justified when the fiber has a much greater heat 
capacity and thermal diffusivity than the surrounding medium. 

3, Oscillating Surface Temperature 

The mathematical analysis is still simple in this case; we 
have just shifted the heat source from a volume to a surface. We set 

A = 0 (A-ll) 

o 

and 


T(r,t) 


r=a 



(A-12) 
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where A is defined only at the fiber surface* The internal and ex- 
a 

ternal temperature fields are decoupled. This case closely approximates 
the temperature oscillations of a fiber in a medium. For more general 
conditions the exact boundary conditions must be used. They are dis- 
cussed next. 

4. Non-Isothermal Surface 

In the situation where both the fiber and surrounding medium 
specific heats and thermal dif f usivities are similar, the isothermal 
surface condition cannot be assumed. The heat flux and the temperature 
at the fiber surface must be equal at all times. The mathematical de- 
finition as to what occurs at the fiber surface is defined by [22,38] 


STjCr.t) 

3r 


= K. 


3T E (r,t) 

9r 


r=a 


r=a 


(A-13) 


and 


T x (a, t) = T E (a,t) . (A-14) 

With these two sets of boundary conditions any simple heat 
conduction problem with a cylinder in a different infinite medium can 
be solved. Next the above solution and boundary conditions are applied 
to the problem at hand. 

C. Internal Heat Transfer 

Heat transfer within the fiber with an oscillating surface 
temperature, Eq. A-12, is described first. In this case only the 
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boundary conditions at the center and surface are needed to solve 
Eq. A-3 (A o =0) . The exact solution is 


T ]; (r,t) 


A 

o 


^ 3/2 ) 

A,(^ 3 ' 2 ) 


e joJt 


+ T 

o 


(A-15) 


Grdber [3l] derived relations similar to Eq. A-15 for convective heat 
transfer to a cylinder. The paper has many theoretical examples and is 
quite complete. 

We set 1^=0 and plot the normalized magnitude of Eq, A-15 

(T -T ) /kt A versus distance for various values of £a in Pig. A-l. If we 
I o a 

define the thermal boundary layer thickness as 6=/2a^7co, we find that 
£=/2 /<5. The internal thermal boundary layer is illustrated by the 
region near the surface which changes its temperature, and, as 
becomes larger, the relative boundary layer thickness becomes smaller, 
indicating that heat transfer is less important when 10* The in- 

ternal thermal boundary layer of the fiber is significant only when 
the frequency is low, the fiber diameter is small, and/or the thermal 
diffusivity is large (£^a< 10). 

D, External Heat Transfer 

Heat transfer outside the fiber with an isothermal surface, 

Eq. A-12, is described next. In this case only the boundary conditions 
at the fiber surface and infinity are needed to solve Eq . (A-3) . The 


exact solution is 
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FIGURE A-1 

NORMALIZED MAGNITUDE OF THE INTERNAL, OSCILLATORY TEMPERATURE 
VERSUS RADIAL DISTANCE FOR A FIBER WITH AN ISOTHERMAL SURFACE. 
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A 

T (r,t) = -2- 
E ju) 




e JU,t + T 


(A-16) 


We set T o =0 and plot the normalized magnitude of Eq . A-16 
( t e - t o )/kA o versus distance, for various values of £ £ a, in Fig. A-2. 
The external thermal boundary layer is illustrated here. As shown in 
part C, when S^a increases the relative boundary layer thickness de- 
creases. The boundary layer thickness is important only when £ a< 1.0. 

E 

E. Coupled Heat Transfer 

This case is not so simple because the internal and external 
fields are coupled. The result is two coupled equations which define 
the total thermal field. 

Inside the fiber the solution is 


Tj(r,t) = 


j 


wE 


e J “ £ + I 


(A-17) 


and outside the fiber the solution is 


T (r,t) = -H 
E u) 


j + 


C I K I J 




3/2 


K(«E r J 1/2 ) 


iu)t . _ 
e + T 


(A-18) 


where 


E = 5. 


^E K o( { E a:| 1 / 2 ) J o( E I a J 3/2 ) + 35 I K I J o( ! I aj 3 / 2 ) K o( 5 E a:lI/2 ) (A ‘ 19) 
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A SINGLE ISOTHERMAL FIBER. 



and the primes denote derivatives with respect to the function argument. 

A plot of these equations versus distance, for several values 
of £ T a and £_a, is shown in Fig. A-3. This plot shows that the iso- 
thermal surface boundary conditions are adequate for our purposes. The 
dotted curve is for a Kevlar fiber in air and a 0.1 kHz temperature 
oscillation. As the frequency increases, the boundary layer becomes 
smaller and the fiber-air system becomes more isothermal, as shown by 
the solid curve. The dashed curve indicates what the temperature oscil- 
lations would be for two materials of similar properties (a Teflon fiber 
in nylon) . 

In this section the oscillatory heat transfer in and around a 
cylinder has been studied. For the purposes of our study it has been 
shown that the fiber material is inconsequential to the formation of 
the thermal boundary layer in air and the isothermal surface boundary 
condition can be used in the analysis in part II of this appendix. 
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II. Oscillatory Heat Transfer Effects in an Air Saturated, 

Fibrous Porous Material 

In the previous analysis the only consideration was the heat 
transfer to and from a single fiber in an infinite medium. In this 
section the conservation of energy equation, which turns out to have 
the form of Eq. A-l, is used to study the heat transfer to and from 
a single fiber when the temperature oscillation is caused by an acoustic 
wave for which X >> a. The results are used to determine the heat 
transfer effects on the phase speed and attenuation in a porous mate- 
rial. The approach in this section is modeled after an analysis by 
Pfreim [63] and restated by Devin [26] on heat transfer effects in 
oscillating gas bubbles in water. 

In the following, the energy equation is derived and solved, 
the sound speed is determined, and then the effects of viscosity and 
heat transfer are simply added and a complete small signal solution 
obtained. 


A. Energy Equation 

The process to be described must obey the first law of 
thermodynamics, the conservation of energy 


dE _ d£ dW 
dt dt + dt 


(A-20) 


where dE/dt is the rate of increase in the internal energy of the gas 
per unit volume, dQ/dt is the rate of heat transferred to the gas per 
unit volume, and dW/dt is the rate of work done on the gas per unit 
volume. Since the fluid considered here is a perfect gas, the rate of 
increase of energy per unit volume is 
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dJ= pC ^ 

dt v dt 


(A-21) 


where C is the specific heat at constant volume. Our analysis is for 
v 

small signals and we can set p = P q , where p Q is the ambient density, 
and use partial derivatives in Eq. A— 21 and the following. The rate 
of heat transferred to the gas per unit volume (in cylindrical coordi- 


nates) is 


Vt ^ 1 3T 

= K ~ H T 

y r 3r 


(A-22) 


The rate of work done on the gas per unit volume is 

aw = p _3p_ 
at ” p at 


(A-23) 


where P is the total pressure. The only overt differences between the 
above equations and the equations in Devin's analysis [26] are that 
Devin was concerned with spherical coordinates and viewed everything on 
a per unit mass basis. Equations A— 21, A— 22, and A— 23 are used in Eq. 
A-20 to yield the energy equation 


„ a t fa 2 T , l 3 t 1 . p afi. 

o v9t ^ 2 r 9 r I p3t 


We need to replace the density term with a pressure. The 
gas is assumed to be ideal and 


P = p(C - C )T 
P v 


where C is the specific heat at constant pressure. If a time 
P 
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derivative of Eq. A-25 is taken and Eq. A-25 substituted into the 
derivative, the result is 


P_ dg_ _ dP 
p dt " dt 




(A-26) 


Substitute Eq . A-26 into Eq . A-24, eliminate like terms, and write the 
equation as 


2 

8 T ,IiT 1 3P 

.2 r 9r a f _ 3t k 9t 

dr I 


(A-27) 


Equation A-27 is a general equation applicable to a small 
signal acoustic wave in the gas whenever X >> a. Equation A-27 is 
equivalent to Eq. A-l, with 

A = -(9P/3t)/K . (A-28) 

The analysis can probably be performed for waves other than sinusoids, 
but, for the analysis here, sinusoids are the only waves considered. 

B. General Solution 

The solution of Eq . A-27 is accomplished by assuming that the 
acoustic wave is a sinusoid such that 

P = P q + pe jaJt , (A-29) 

where P q is the ambient pressure, p is the magnitude of the acoustic 
pressure, and X >> a. Since X >> a and this is a boundary layer problem, 
there is no pressure gradient in the vicinity of the fiber; in other 
words the acoustic pressure p is assumed to have no dependence on the 
radial distance r. Equation A-27 becomes 
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3r T 


(A-30) 


If we assume that the fiber temperature remains constant, then Eq . A-16 


is used to obtain the solution of Eq . A-30. The solution is 


T(r > C) " T o + fc 
o p 


K Q (CrjS 


K o (Saj^) 


This solution is used in the next section to calculate the heat 
transfer effects on the sound speed in a fibrous porous material. 

C. Heat Transfer Effects on the Phase Speed 

Now that the radial dependence of the temperature is known, 
the small signal phase speed may be calculated. Since we have losses 
caused by heat transfer effects, the phase speed is neither adiabatic 
nor isothermal. The phase speed is written as 


2 do 
C T dp 


(A-32) 


The limitation that the compressions be isentropic has been removed 
because the compressions are no longer lossless. From Eq . A-32 the 
phase speed can only be calculated if a relation between the pressure 
and density can be found. In order to do this the heat transfer 
effects are integrated from the fiber surface out to some reference 
distance. The specific reference distance is determined later in this 
appendix. The volume integration is performed in cylindrical coordi- 
nates. 

In cylindrical coordinates the specific volume per unit mass 
and per unit length of any shell is cylindrical 


! 


1 


] 


1 
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— = 2-rrrdr 
p o 


and infinitesimal changes in pressure and temperature are written, 


respectively, as 


dP = pe° 


dT = 0 (r)e J 


The relation of the thermodynamic states between any two shells is 


p p T ’ 

H o o 

The above equations are now combined to determine a differen- 
tial form of Eq. A-31. Differentiate Eq . A-36 and substitute Eqs. A-33 
through A-35 and Eq. A-31 into the differential of Eq. A-36 to obtain 


d( I) = _ ilrpe_ 
P YP„ 


K (SrjS 

1 + (y-1) 1 - — xr dr ■■ 

K (Sajp) 


Integration of Eq . A-37 is accomplished by noting that [5 9 .] 


xK (xj^)dx = - -t-K (xj' 5 ) 
o . ^2 1 


= - xj (kei^x + j kei_^x) 
= x(kei T x + j kei T x) 


The upper limit of integration in r is taken as arbitrary and the exact 


value is determined later. After integrating over the volume around 
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the fiber and some algebra the resulting equation is 


F - P 
o 

YP„ 


2(y- 1) j ^[ri^CCrj’ 5 ) - aK^Saj' 5 )] 

C(r 2 -a 2 ) K o Uaj h ) 


1- 


(A-39) 


The last step is to take the derivative of this equation and solve for 
dP/dp. This result yields 


c 2 _i!o 

Ct p o 


1- 


2(y~ 1) j S5 [rK 1 (5rj ,s ) - aK^a;^)] 
5 (r 2 - a 2 )K Q (Caj^) 


-1 


(A-40) 


Although desirable, the transform to the time domain could not be 
accomplished. Thus, this equation, in its present form, is only good 
for single frequencies. 

It is important to determine whether the function has the 
proper values in the limits of fiber packing. If the fibers are spaced 
infinitely far apart, then r+°° and 


2 ^ P o 


= 


T p. 


(A-41) 


which is the square of the adiabatic sound speed. 

In the limit in which the fibers are packed as close together 
as possible, r+a, and the equation becomes 


? P 
l o 


(A-42) 


which is the square of the isothermal sound speed. 

In a porous material the fibers are packed in a spacing well 
between these two extremes. The relative fiber packing is used to 
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determine the upper limit of integration. Assume that the fibers are 
evenly spaced and parallel, as shown in Fig. A-4. The dotted lines 
indicate the radius at which the thermal boundary layers intersect. 

The small spaces labeled M indicate the area missed by the integration 
(19%), but this area has little effect on the calculated phase speed as 
Eq. A-45 is a slowly varying function. 

For two cases, P = 0.918, the magnitude of the sound speed is 
plotted versus frequency in Fig. A-5. Since the largest phase angle of 
Eq. A-45 is less than 2°, the phase angle is ignored in further cal- 
culations. The magnitude change of the sound speed is found to be a 
very slowly varying function of both porosity and frequency. Increase 
of the thermal diffusivity will cause the phase speed to approach iso- 
thermal. This effect is the same as lowering the frequency, since 
£ = AT/a^ . For a constant porosity, reduction of the fiber size causes 
the fibers to be more tightly packed, as shown in Fig. A-4 (r*a) , 
Increased packing density causes a lower phase speed. 

It is now important to compare the heat transfer effects to 
the viscous effects in the same porous material. This is the topic of 
the next section. 

D. Thermoviscous Effects on the Phase Speed 

In this section the heat transfer and viscosity effects in a 
porous material are combined. Since only low- intensity sound is being 
considered in this appendix, a simple replacement of b^ by Eq. A-40 in 
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FIGURE A4 

SKETCH OF THE FIBER ARRANGEMENT IN AN IDEAL POROUS MATERIAL. 
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FIGURE A -5 

MAGNITUDE OF THE SOUND PSEED VERSUS FREQUENCY FOR TWO 
POROSITIES OF A FIBROUS POROUS MATERIAL. 




Eq. 1-16 is assumed to be valid. The result of the above substitution 
plotted versus frequency in Fig. A-6. 

In Fig. A-6 the thermoviscous phase speed predictions and 
data for Kevlar 29 (p = 0.980) are plotted versus frequency. This 
combination works fairly well for the full frequency range. 




FIGURE A -6 

MEASURED AND PREDICTED THERMO-VISCOUS PHASE SPEED VERSUS 
FREQUENCY FOR BATTED KEVLAR 29 (P = 0.980). 






III. Summary and Conclusions 


In this appendix the oscillatory heat transfer effects in fi- 
brous porous materials have been determined. The results have been com- 
bined with viscous effects and used to predict phase speed in a fibrous 
porous material. In the analysis, the fiber type was found to have 
minimal effect on the heat transfer boundary layer in air surrounding 
the fibers and an isothermal boundary condition is used in predicting 
the soundspeed. 

As anticipated, the heat transfer effects on the sound speed 
are bounded by the isothermal and the adiabatic sound speeds. These 
limits encompass a change of about 18% from the adiabatic sound speed. 

Because of the complication of using Eq. A-45 in the body of 
this dissertation, the limiting cases of the isothermal and adiabatic 
sound speeds are used in the low- intensity sound predictions. When 
intense sounds propagate through a porous material the acoustic particle 
velocity and displacement become large. The large particle displace- 
ment causes the heat transfer to be convective, rather than conductive, 
and causes acoustic streaming which mixes the fluid. Both of these 
effects shift the curves in Fig. A-5 so that, for a given frequency u), 
the fluid compressions are more isothermal and the sound speed is 
slower. In addition, the isothermal sound speed is used in the high- 
intensity sound predictions because of the increased heat transfer and 
gas mixing caused by the high particle velocities. 
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APPENDIX B 


DRAWINGS 

Construction drawings of two devices built for the experiments 
in this study are presented in this appendix. Figures B-l through B-3 
show the construction drawings for the standing wave tube. Figure B-4 
shows the construction drawing for the traveling wave tube. 

In Fig. B-l -A the impedance tube, the sample holder, and the 
aluminum plug are drawn. In Fig. B-l-B the driver/probe tube adaptor 
for the impedance tube is drawn. In Fig. B-2 the impedance tube 
support/carriage track is drawn. In Fig. B-3 the microphone carriage 
and microphone holder are drawn. 

In Fig. B-4 the traveling wave tube and its accessories are 

drawn. 
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STANDING WAVE TUBE. 
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APPENDIX C 


COMPUTER PROGRAMS 

Two of the programs used in this study are listed here. A 
short description is given for each program. 

PERT4PD 

This interactive program is used to calculate the propagation 
of an initial wave consisting of a fundamental and a second harmonic 
component of arbitrary amplitudes and relative phase in a porous 
material. The program calculates and plots the sound levels of the 
four harmonic components with distance (0-25 cm) . The information 
needed to use the program is the frequency f, linear flow resistivity 
a/415, nonlinear flow resistivity q/415, fundamental component sound 
level, porosity P, second harmonic component sound level, and relative 
phase <j>. 

PHSPED 

This interactive program is used to calculate the heat 
transfer effects on the phase speed of sound in a porous material. The 
program calculates and plots the magnitude and phase angle of the phase 
speed versus frequency. The information needed to use the program is 
the thermal conductivity, density, and specific heat of the fibrous 
material, the fiber radius, and the half spacing between the fibers. 
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PROGRAM PERT4P0 (INPUT t OUTPUT *DATPTSt PLOT, TAPE5=PL0T) 

c 

C HLK (PAD) 6-1-B2 

C THIS PROGRAM IS USED TO CALCULATE THE AMPLITUDES OF THE FIRST 

C FOUR HARMONIC COMPONENTS, WTTH RESPECT To DISTANCE* FOR AN 

c acoustic wave traveling in a bulk pohoijs material. 

C THE BOUNDARY condition CONSISTS OF AN ARRTTRARY COMBINATION 

c OF THF FIRST and SECOND HAHmONICS. DATA FOR THE FIRST FOUR HftRMONlC 

C COMPONENTS ARE PLOTTED. IhE IMAGINARY PART OF THE PROPAGATED WAVE 

C AMPLITUDES ARE PLOTTEC. 

C 

INtEGFR LABX ( I > fLABYl (1 ) *LarY2 ( 1 ) ,MAT (2) 

REAL PlPO (1601 )»X( 1601) »P2PO(1601) 

REAL PlP02(l6ol) ,P3P0 ( 1601 ) .P4POU601) 

REAL XD<21) *FnD<21> , 2ND <21 > .TPO<21> ,FRTH(?1) 

COMPLEX 01 ,Q2,03,0A,U2l *U2^,U23»U2A.V1 ,V? 

COMPLFX GAMMA 1 ,GAMMa2, J,GAMa 1C* GAMA2C,GAmmA3, gamma* 

C 

C SPECIFICATION OF CONSTANTS FOp PrOGrAM, 

C 

AXLEN=B,5 

AYlEN*5.0 

J*CMPlX(0.,1.) 

C=2V2. 

RH00=1.2l 

PI*0.0 

LDAT*6LDATPTS 

c 

5 CONTINUE 

c 

C DATA entry 

c 

PrINT**MNPUT NO. POllNTS/lNCH TO PLOT IlfiO MAX.l, IOPT ** 

reau»,npinch»iopt 

c 

c program input termination. 

c 

IF (NPINCH.EO.O) stop 

c 

C DATA w = OF DATA LlNEs» p LOT *TITLE*fPOBOSlTY: DECIMAL* SlG M A» LINEAR 

C COEFFICIENT OF RESISTIVITY PER RHOO»C, ETAS NONLINEAR COEFFICIENT OF 

C RESISTIVITY Per RHO0*C* frequency: hertz, B: relative 
c amplitude of the second harmonic component with respect 
c to the fundamental* psi: Rti ative phase between the 
C FUNDAMENTAL AND SECOND "HARMONIC COMPONENTS. 

c 

R£AD(lOAT»*)NO*MAT(1) »MAT(2) *P.F 
IF(ND.EQ.O) GO TO 500 

c 

OMEGAs6.2B318S*F 

c 

print**f$Igma, eta* p$i+f* 
reau*. sigma* eta* psi 
c 

DO 150 JD=l,NO 

READfLDAT**) XD( JD) *FiND( JD) ,SN'0( JD) *TRO( jn) ,FRTH( JD) 

IF(JO.GT.l) go TO 17 
EPs*lO**< (FND(l)-lR 1 . 2 )/ 20 .) 

B=10««( (SND(l)-FND(l) )/20.) 

7ERO=FNO ( 1 ) 

17 FNoUn)=-ZERO*FND(JD> 
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5Nq< J n)=-ZEP0.SND<JD> 

TRoUn)=-ZEPO*TRO(JO) 

FRTH< JD> *«ZERo*FRTH ( JO 
IF(SNr><JD) .LT.- 50 . 0 ) SND(JO)* 10.0 
IF <TRn< JO>.LT.- 50 . 0 ) TRD( JD)= 10.0 
IF(FRTH(JD) .LT.- 50 . 0 ) FRTH ( jD) = 1 0 . 0 
150 CONTINUE 

c 

SlGMAl=415,«SIGMA/(RHCO*OMtr,A) 

ETaIsAIS.^ETA^C/ (RHOO*OMEGA) 

gamma 1 =CSQRT ( 1 ,-J« 5 lGN Al ) 

GAMMA?=CS 0 RT ( 1 ,-J*SIGKA 1 / 2 .) 

GAMA 1 C»CSQRT< 1 .+J*SIGPA 1 > 

GAMA 2 C=CSQRT (l.*J*SlGMAl/ 2 .> 

GAMMA 1 =CSQRT ( 1 .-J*SlGPAl/ 3 . ) 

GAMMA 4 =CSQRT <l.-J«SlGNAl/ 4 .) 

NPsNPTNCH*AXLEN 
DElTAX= 25 . 0 * 0 mEGA/ <NP« 100 .*C> 

SVaL=- 10 .*AYLEN 
PSl* 0 . 0 l 74 S«PSI 

computation of the first four harmonic component ampituDEs* 

WITH RESPECT to DISTANCE* FnR AN ACOUSTIC WAVE TRAVELIING 
IN A RULk'pORoUS MATERIAL. 

00 10 JL= 1 ,NP 
X(JL)=OElTaX*JL 

IF noPT.EQ.l) GO TO IS v , , , , 

IF (IOPT.EQ. 0 ) P 1 PO(JL)= 20 .*ALOG 10 (CABS(CEXP(.J*GAMMA 1 *X(JL) ) 1 ) 
IF dnPT.EO.O .AND. PlPO ( JL) .LT.SVAL* PlPO<JU=SVAL 
15 V 1 =B* ( 2 .*GAMMa 2 -GAMA 1 C) 

V2=Z.*ETA1-J*(GAMMA1 # «2) * (GAMMA24GAMA1C) 

V 1 =V 1 »V 2 

V 2 =o.*GAMA 1 C*GAMMA 2 » U ,-GAMa 1 C*GaMMA 2 ) 

01 =V 1 /V 2 

VI =J*fTA 1 -J* 2 .*( SIGMA 1 *J* 2 .)*GAMMA 1 
02 =V 1 / ( 2 .*SIGma 1 *GAMMa 1 ) 

V 1 =ETa 1 ♦J»GAMma 1 * ( 3 ,hl*SIGma 1 ) 

V 2 = 4 ,*GAMMa 2 *GAMMA 1 * ( 1 ,-GAMmA**GAMMA 2 ) 

Q 3 =B* (2.*GamMa2*GAMmA 1 ) *Vl/v 2 
Vl= (B*» 2 ) f < 2 .*<GAMMa 2 <, * 2 >*StGPA 1 > 

V 2 =J«fTA 1 -J* 4 .*GAMMA 2 « (SIGma 1 *J* 4 . ) 

Q4=V1*V2 

vucexp <j*( -Gamma i*x C v l> *pst ) ) 

V2=CEXP < J* (- <2.*GAMMA2-GAMaiC) *X < JL) *PSI )) 

U 2 l* ( V 2 -V 1 ) 

V 1 =CEXP( J* 2 .*(-GAMMA 2 *X(JL) ) ) 

V2=CEXP(J*2.*(-GAMMA1«X(JL) ) ) 

U22*(V2-V1) 

VlaCExP < J* (3.* (-GAMMA3*X ( JL) ) ♦PSI ) ) 

V 2 =CE*P < J* (- ( 2 .*GAMmA 2 »GAMma 1 ) *X ( JL) *PSI ) ) 

U 23 *(V 2 -V 1 ) 

V1=CEXP(J*(4.*(-GAMMA4*X( JL) > *2**PSI) ) 

V2*CExP ( J # ( 4 .* ( -GAMMA 2 *X ( JL D +2 • *PSI > ) 

U 24 =(V 2 -V 1 ) 

V 1 =CEXP(-J*GAMMA 1 *X(JL) ) 

V1=V1.EPS*Q1*U21 

V2=B*CEXP(-J»(2.*GAMMA2*X( J| ) ♦PSD ) 

P 1 P 02 ( JL) * 20 .«ALOG 10 (CABS ( < V 1 ) ) ) 

P 2 PO<JL>= 2 o.*ALOG 10 (CABS( ( V?+EpS* 02 *U 22 > ) ) 

P 3 PO(JL)= 2 o.*ALOG 10 (CABSf <EPS*Q 3 *U 23 > ) ) 

P 4 P 0 < JL > = 2 o.*AL 0610 (CABS ( (EPS*Q 4 *U 24 > 1 ) 
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IF IP1P02 ( JL> .LT.SVAL) PlPo?(JL)aSVAL 
IF (P?P0 ( JL> .LT.SVAL) P2P0<JL)=SVAL 
IF (P3P0 (JL) .LT.SVAL) P3P0 ( JL ) =SVAL 
IF (P4P0 (JL> .LT.SVAL) P4P0( JL)=SVAL 
lo CON I InuE 

c 

DO 30 JN=i,KP 
X ( JN) =X ( JN) *100.*C/0MEGA 
30 CON I INUE 

c 

c**** PLOT THE VALUES. 

C 

L Apx n ) -*x . cm* ] 

LAryi (1)=*l0G(U1/U0>* 1 

LARY2(l)s*LOG(UI/U0)* 

xopg=i .0 
yopg=?.o 

DXa25.0/AXLEN 

DYalO. 


branching for plotting choices. 

PLOI Ml/uo 

IF (IOPT.E0.1) GOTO 20 } 

call PLTLFN(4LPL0T) 

call PLT 0 lM(ll . 0 * 8 . 5 * 1 * 4 . 0 ) 

CALL pltorg (XORG.YORG) 

DRAW AXES 

call OLTAXIS <0 . 0 * 0 . 0 * AXLEN.o.O * 0 . 0 »?s. 0 * 1 . 0 *LABX.- 6 * 5 *-* 1 *-.l ) 

CALL PLTAXIS < 0 . 0 * 0 . 0 * AYLEN.R 0 . *-S 0 . * 0 . 0 *?. 5 »LABY 1 , 10 * 5 *-. 1 *-.l ) 
call PLTAXIS ( 0 . 0 .AYLEN .AXLtN* 0 . 0 * 0 . 0 * 25 . 0 * 1 . 0 *LA 8 X.- 6 * 0 * . 1 *.l ) 

CALL PLTAXIS (AXLEN.O.O.AYLtN. 90. 0.-5V..0, 0*2. 5.LABY1, 10* 0..1..1) 

PLOT THE DATA. } 

i 

CALL PLTDAtA (X.P1P0.NP, 0,0.0. ,DX,-50.0.Dy,.08*2) 
call PLTLInE <4.25»-1.2»"0.14) 

WRITE (5* 100) P.F, SIGMA 
CALL PLTENDUl .0*8.5) 

20 CONTINUE 


plot uiuo. 

CALL PLTlFn(4lPL0T) 

CALL PLTDIM(11.0*8.5,1,4.0) 
call PLTORG (XORG.YORG) 

C DRAW AXES 

CALL PL TAXJS <0.0*0.0*AXLEN. 0.0*0.0*25.0*1 . 0 *LABX * -6 *5*-* 1 * -. 1 ) 
CALL PLTAXIS (0.0.0.0* AyLEN.qO, ,-50. ,0.0*?.5,LABY2, 10,5*-. 1 1) 
CALL PLTAXIS (O.O.AYlEK, AXL tN, 0.0,0. 0*25. 0.1. 0,LABX, -6,0* .1 . .1 ) 
CALL PLTAXIS ( AXLEN* 0*0* A YLtN* 90.0,-50. *0.0*2. 5.LABY2* 10*0. .1 * .1 ) 


l 
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CALL PLTOAJA (X,PlP02»INP*0tO,0,nX,-50.0.DY,.08,2) 

CALL dLTDATA (X.P?P0.NP.0.0.O.DX<-5O.O.DY,.O8.9) 

CALL PLTOATA <X,P3P0.NP.O.O.n»DX.-5O.0.DY,.O8.4> 

CALL PLTOATA <X»P4P0.NP.0.0.O*DX*-50.0»DY,.08.8> 

CALL PLTOATA <XD.FND.NC.-1»?.0.,DX.-50.0»OY*0.16»3> 
CALL PLTOATA <XDtSND*NC*-l*5.0..DXt-50.0*DY*0.l6t3) 
IF(TRn(l) .GT.10.0) GO TO 6b 

CALL PLTOATA (XD.TRD.NC.-l * 9 . 0 . . Dx .-50 .0 *DY .0 . 16 *3 > 

65 IF(FRTH(1) .GT.0.0) GO TO 7b 

CALL PLTOATA (XD»FRTH»iNO«-1 «1 1 *0»»DX,-50»0»DY.0*16»3> 
75 CALL PLTLINE <4 .25 * -0 »9 . -0 . 1 * > 

WRTTE(5*200> p.F.ZERO 
CALL pLTlINE(4.25,-1.3.-0.14) 

WRt TE (5*300' SIGMA. ETA 
call plTLINE (4. 25. -1.8. -0.14) 

WRj TE (5*400 ) RtPSI 
CALL PLTENDUl.O. 8 . 5 ) 

plot labeling for inRlt parameters. 

100 FOrMAt (AP = A.F9.3.A F = A.F9.3.F SIGMa s *,F9.3> 
200 FORMATIAP = A.F9.3,* F * F.F9.3,* SPL = A.F9.3) 

300 FORMAT (ASIGMA = *,F9.3t* ETA r A.F9.3) 

*00 FOPMAT<*B = A.F9.3.* PSI s A.F9.3) 

GOTO 5 

C 

500 PR I NT* .MAT ( 1 ) .MAT (2) .* SAVE PLOT S* 

c 

ENo 
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PROGRAM PHSPEf) (INPUT, 0 UTPut»PL0T»0UT,TApe 2=0UT,TAPE5=PL0T> 
COMPLEX C2S.C2S1 »C2sC,ViVl , v2,v 3 
INf tGFR LABX (?) ,LABy (2) «LAby2 (2> 

REAL T2SM( 1701 ) ,L0GF (1701 ) , PHASE ( 1701 ) ,F ,REP, AEP ,0LF ,KE«A , KE I A » 
* KERDA,KElOA,KERDP,KEIUR,KI ,PI »CPI 

INITIALI7E VARIABLES 
GAmmAsI .402 
ALPHAF=2.216E-5 

REp1=SQRT(2.«3.1415927/ALPhaE) 

AEP1=REP1 

S PRINT*, /ENTER KI, PI, CPI/, 

REau*,KI ,PI ,CPI 

IF(PI.EQ.O.O>STOP 

PRINT*, /ENTER FIBER RADIUS, RADIUS, NO PLOTTER STEPS/, 
REau*,A*R»N 

REP=RFP1*R 

aep»afpi*a 

0LF=2./N 

ALP"AT=KI/PI/cPI 

AlpaSORT (2,«3. 14 15927/ ALPHA I ) *A 
VSR=SORT(31.0711/(KI«PI*CP1) ) 

M=n*1 

COMPUTE PHASE SPEED. 

DO 10 J=1 «M 

LOGE ( J ) =2. ♦ ( J-l ) *DLF 

F=lO.** (LOGF ( J> ) 

RE=HEP*SQRT (F) 

AE=AEP*SQRT (F) 

AI=AIP*SQRT (F) 

call MMKEL0 (AE*BER*BEI,KERa,KETA,IER) 
call MMKEL0 (Ae,BER,REI,KERua,KEIDA,IER) 
call MMKELO (RE,BER»BEI ,KERuR,KEIDR,IER) 
call MMKELO < a 1 ,8ERA .BE I A , XkfR , XKE I ,IER) 

CALL MMKELO (AI ,BERDA»BEIDA,xKeR,XKEI«IER) 

XlrHE*KElDR-AE*KEIDA 
X2=AE*KER0A-RE»KERDR 
C2S1=CMPLX (XI ,X2) 

C2SU=CMPLX (KERA.KEIA) 

V 1 =CMPLX (KEROA.KEIDA) 

V2=LMPLX (BERA.BEIA) 

V 3 rCMPL X (BERDA.BEIDA) 

V=VSR*Vl*V2/V3-C2SO 

C2S=343./CS0RT ( 1.-2* (GAMMA-i . >*C 2S1/ ( (RE*RE-AE*AE > *V ) ) 
C2SM(j)=CABS(C2S) 

PHASE (J> *ATANp (AIMAG(C2 S) ,HEAl (C2S) ) *180 . /3 . 1415927 
10 CONflNUE 

Print out values. 

WRI I E ( 2, 1 50 ) 

DO 20 J=1 ,M 
F=10.**L0GF ( J) 

WRITE(2»*)F.C2SM(J) .PEASE (J) 

20 COnIINUE 
150 FORMAT(AlA) 
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c plot magnitude and phase vs. Frequency, 

DATA LA8X//FRFQUENCY,/,/ HERTZ*/, LABY//NORM PHASE/*/ SPEED// 
DATA LA8Y2//PHASE ANGL/./E// 

XORlis 1 .0 

yorg=?.o 
A XleN* 8 • 5 
AYLEN=5.0 
DX=2 ,/aXlEN 

DYrAOO./AYLEN 

DYj=2o./AYLEN 

c 

CALL PLTLEN <4 lPL 0T> 

call PLTDIM<U. 0*8. 5*1*4. 0) 

CALL pltorg <X 0RG*Y0RG> 

C draw axes. 

call LOGAXIS < 0. 0*0.0 * AXLEN, o.0«100.,1000n. *LARX,-] 6*-. 1 ) 

CALL PLTaxis ( 0.0*0.0*AYLEN,g0.0*0.0.A00.,50.»LABY.16.2»-.l «-.l ) 
CALL LOGAXIS <0.0 * A YlEK, AXLEN , 0.0,1 00. *10 ooO.,LABX, 0*-,1 ) 

CALL PLTAXIS (AXLEN *0.0 *AYLEN*90«0 *0.0*1 .0*.! *LABY* 0*0 ».l *.l) 

c 

call PLTDAtA (L 0GF.C29F *M*0*o*2, *DX*0*DY*n.) 

CALL PLTlINE( 2.0*-1.2,.14) 

WRI IE ( 5 * 1 00 > A,R 
CALL PLTEND (11,0*8.5) 

plot phasf angle. 

CALL PLTlFn(AlPLOT) 
call PLTDlMdl. 0*8. 5*1*4. 0 ) 
call PLTORG(XORG.YORG) 

Draw axes, 

call LOGAXIS ( 0.0 *0.0 * AXLEN.o.O* lo. * 10000., LABX*-16*-.l) 
call PLTAXIS ( 0.0*0.0*AYLEN,90.0*-10,*10.*1.*LABY2*11 ,2*-.l ,-.l ) 
CALL LOGAXIS (O.O.AYlEF, AXLEN* 0.0, 10. *10000.. LABX.O*-. 1 ) 

CALL PLTAXIS ( AXLEN. O.o . AYLEN *90 • 0 *-l 0 . ,ln. *1. ,LABY2*0*0* .1 . . 1 ) 

call PLTDATA (LOGF,PHAsE*M*O,0, 1 . ,Ox,-10.,DY2*0.0 ) 
call PLTEND ( 11.0,8.5) 


GOTO 5 

c 

loo FORMAT <*A-/,E9.3,10X*/RT*.t9.3) 

END 


1 


1 
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